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Abstract 



Wc compute the distribution of the purity for random density ma- 
trices (i.e. random mixed states) in a large quantum system, distributed 
according to the Bures measure. The full distribution of the purity is com- 
H puted using a mapping to random matrix theory and then a Coulomb gas 

method. We find three regimes that correspond to two phase transitions 
(~j in the associated Coulomb gas. The first transition is characterized by an 

Q explosion of the third derivative on the left of the transition point. The 

O second transition is of first order, it is characterized by the detachement of 

a single charge of the Coulomb gas. A key remark in this paper is that the 
I random Bures states are closely related to the 0(n) model forn = 1. This 

^ actually led us to study "generalized Bures states" by keeping n general 

instead of specializing to n = 1. 
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1 Introduction 

1.1 Why random states ? 

In the last decades, there has been an increasing interest in random states in 
quantum physics. Taking a state at random corresponds in some sense to assume 
minimal prior knowledge about the quantum system [TH]. Random states are 
thus very general, they can be seen as typical states. They play the role of 
reference states to which the evolution in time of a given physical state can be 
compared. The idea behind random states is close to the one of Wigner when 
he introduced random matrices to describe the Hamiltonian of a large nucleus 
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-a very complex system with many unknown details. In a similar way, in a large 
quantum system when the Hamiltonian is not known precisely one can model 
the state of the system by a random state. 

Let us consider a quantum system described by a Hilbert space % of dimen- 
sion TV. A pure state is a vector G T-L of unit norm, it describes the state of 
the quantum system. There is a natural way of choosing pure states at random: 
it is to pick them according to the uniform distribution. The uniform distribu- 
tion -also called Haar measure- is indeed the unique probability measure that is 
invariant under a change of basis (unitary transformation) . Random pure states 
have been studied extensively [nilHIIlllSSlllllSllllSl- In particular there has 
been a lot of studies about entanglement properties of such random states in a 
bipartite system [HI [331 [231 US HOI 132] j i-e. when the system is made of two 
subsystems V. — Ha ^Hb- 

The notion of quantum state can be generalized by introducing some statis- 
tical uncertainty. A generalized quantum state or "mixed state" is described by 
a density matrix cr, that is a iV x A'^ matrix on % with the following properties: 
-it is Hermitian, i.e. — a 

-it is semi-definite positive, i.e. (u|(t|u) > for all \v) e 
-it has trace 1, i.e. Trcr = 1. 

The matrix a can thus be diagonalized in an orthonormal basis {|wi)}i<i<Ar, it 
has A'^ nonnegative eigenvalues that sum up to one. The spectral decompo- 
sition of (T is thus of the form: 



The eigenvalue can be interpreted as the probability to find the system in pure 
state \vi). cr is a statistical superposition of pure states. If only one eigenvalue, 
say Aio, is non-zero, then Ai^ = 1 and a = \vig){vi^^ \ is the projector on \vig), a is 
thus simply the density matrix representation of pure state lui^) -we will call a 
a pure state in this case. If several A^ are non-zero then there is an uncertainty 
about the state of the quantum system. In general, such a mixed state can 
describe a system that is not isolated but in contact with an environment. The 
interaction with the environment, more precisely the entanglement between the 
system and the environment, is responsible for a loss of information for the 
observer that can do measurements only on the system. The more mixed the 
state of the system is, the more entangled the system is with its environment. 

One can now be interested in taking at random mixed states. There is no 
single natural probability measure for mixed states (unlike pure states), as the 
unitary invariance only implies that the measure depends on the eigenvalue 
spectrum alone (thus not on the eigenbasis of the density matrix). Several 
probability measures have been proposed for mixed states and motivated on 
physical grounds. One possibility is to construct a measure from a unitarily- 
invariant distance between mixed states. Such a distance induces a measure on 
the space of density matrices (volume element associated with the distance). 
By construction, this measure is invariant under unitary transformations, it is 



N 



N 




(1) 
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thus of the form P(Ai, . . . , AAr)dAi • • • cIAat x df where dv is the Haar measure 
on the group of unitary matrices U{N). The measure is thus characterized by 
the marginal distribution of the eigenvalues V{Xi, . . . , Aat). 

A natural choice of distance is the Hilbert-Schmidt distance [T5] defined by: 

DHs{p,<Tf = TY{{p^af} (2) 

When both density matrices are pure states, i.e. p = and a ~ \(j>){<j)\, then 

p^ — p, — a and the distance reduces to D^s {Pi'^)^ = 2 — 2 The 
Hilbert-Schmidt distance induces a probability measure for mixed states. We 
thus obtain random states that we will call "random Hilbert-Schmidt states". 
The eigenvalues of the density matrix of such random mixed states are dis- 
tributed according to the law: 

Vhs{M.-...\n) = AM5\y^^^-n \{{>^J~>^kf (3) 



where = V{N^)/ Jlfjo' r(iV - j)T{N - j + 1) 



This distribution can ac- 



tually also be obtained in another natural way. If we consider an auxiliary 
space He oi same dimension N (the "environment" ) and choose a pure state 
1-0) at random (uniform Haar measure) in the bipartite space H (g) He, the 
eigenvalues of the reduced density matrix of our system a = Tr^; will 
be distributed exactly according to the law in Eqn ([s]). Tr^ is the partial trace 
over the environment The mapping to random matrix theory for Hilbert- 
Schmidt states can be guessed from the presence of the Vandermonde determi- 
nant det(A^~"^) = Y\i^j{Xj — A,:) in Vhs- The distribution 7'hs(Ai, Ajv) is 
actually the same as the distribution of the eigenvalues of a so-called Wishart 
matrix but with additional constraint ^ • A^ = 1 il3D] . 

A second physically relevant choice of distance is the Bures distance, and 
this is the topic of this article. It is defined by 



Db {a,pf ^2-2^ F{<j,p) 



(4) 



where F {a, p) is the fidelity: 

^/Fi^^TT[[^a^f^] (5) 

The distance Db{<t, p) = Db{p, (j) is nonnegative, it is zero if and only \i p = a. 
When both density matrices are pure states, i.e. p — and a = \(t>){(j)\ and 

thus — p, it reduces to {a, p) — \ and Db{<J,p)^ = 2 — 2 |(V;|0)|.It 

is the topic of this article. The Bures metric is a metric in the mathematical 
sense for density matrices. This choice has been motivated on both measurement 
and statistical grounds [51 [3 [TOl HI] . These random mixed states will be called 



If p is a density matrix on ® "He, then cr = Tr^p = X]^i("^fc|p|"'fc) where {\m^.)}^. 
is a basis of 'He, and cr is a density matrix on H. 
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"random Bures states" . The distribution of the eigenvalues of the density matrix 
for random Bures states is given by [351 HI] • 



N 



n,A 



1/2 



n 

3<k 



Xj + Xk 



(6) 



where the delta function enforces the condition Xi — Ti a = 1 and where 



Cn is a normalisation constant: Cn = 2^'-^ r{N^/2)/ tt^/^ j-j^^^ Y{j + 1) 

Note that the same Vandermonde determinant as in Eqn. ^ appears in ([6) 
but there is an additionnal term in the denominator of the form ni<j (Aj + Ai) 
Osipov, Sommers and Zyczkowski explain in |31j how to generate numerically 
Bures density matrices. 



1.2 Statistical properties of random mixed states: purity 
distribution and entanglement properties 

One can ask several questions about statistical properties of random mixed 
states. The most natural one is to ask how mixed (or pure) such states are. For 
a random mixed state described by a density matrix cr, the purity defined as 
E2 = Tr (T^ = J2i Xi provides an indicator to answer this question. The higher 
the purity, the purer the state. As we have seen above, a purer state corresponds 
to a state that is less entangled with its environment. 

For Hilbert-Schmidt random mixed states, the distribution of the purity has 
been computed recently by one of the authors in collaboration with Majumdar 
and Vergassola for a large system [291 HO] , it features two interesting phase tran- 
sitions. For random Bures states, the average of the purity has been computed 
by Sommers and and Zyczkowski |36| , and the moments were computed by the 
same authors and Osipov [3T] . However the distribution of the purity for Bures 
states was up to now not known analytically, and the goal of this paper is to 
compute it for large systems CH ~ C^, N 00). We find phase transitions that 
are very similar to the ones appearing in the Hilbert-Schmidt case, but with im- 
portant changes in the analytical expressions. The order of the first transition is 
different for Bures states and Hilbert-Schmidt states, but the second transition 
is identical in the two cases: it is a first order transition that corresponds to 
a sudden jump of the maximal eigenvalue. Another very interesting question 
is to consider random states on a bipartite (or even multipartite) system, ie 
H — Ha ®T~Lb (of dimension N = N^Ng), and to study entanglement between 
the subsystems A and B. A mixed state described by a density matrix a on H 
is said to be separable if it can be written as a convex sum of tensor products, 
ie as cr = J2k Pk^^A ® '^b with pt > and J2k Pk — ^- Otherwise it is entangled. 
For mixed states, unlike pure states, one cannot use Von Neumann or Renyi 
entropies as a measure of entanglement, but several measures of entanglement 
have been introduced and studied (such as the entanglement of formation or a 
distance to separable states) [571 133 but they are all quite difficult to compute. 
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especially for random states. An easier question is to ask what is the probabil- 
ity that a random mixed state is separable, or what is the probability that its 
partial transpose is positive (which is implied by separability) [34] . Those prob- 
abilities have been studied by Aubrun et Szarek ^ |37] for the Hilbert-Schmidt 
case, and by Ye |31j for the Bures measure. To summarize the conclusion of 
[Hj, random Bures states are almost surely non separable when iV — >• oo, even 
when conditioned to have a positive partial transform. 

1.3 Relation with the 0{n) matrix model 

A key remark in this article is that the distribution in Eqn ([6| is also the same as 
the distribution of the eigenvalues for the 0{n) matrix model, for the special case 
n = 1, except for the additional constraint J^i'^i — 1- And, at least for large 
N results, this constraint can be handled with help of Lagrange multipliers. 
The 0{n) model was originally introduced by Rostov [52] in a very different 
context, namely quantum gravity and enumeration of random discrete surfaces 
with self- avoiding loops (see |16| for a review of quantum gravity techniques). 
The partition function of the 0{n) model is: 



where the matrices M and Ai are Hermitian N x N matrices and V{M) is for 
example a polynomial in M. The Gaussian integral over matrices Ai can be 
performed. Then M can be parametrized by a unitary transformation and its 
eigenvalues and the integral over the unitary transformation can also be done. 
Finally, Z becomes an integral over the eigenvalues A^: 



In this form, the model can be extended to real values of n (not necessarily 
integer). For the special case of a linear potential V{X) = A and for n = 1 (case 
of a two- matrix model), if we add the constraint X]j = 1, we recover the 
joint distribution of Bures eigenvalues in Eqn ([6|. We can also remark that this 
provides a realization of the Bures measure as the probability measure induced 
on eigenvalues in a two Hermitian matrices model. 

In this paper, we actually study "generalized Bures states" , for which the 
eigenvalues of the density matrix a is distributed according to the probability 
density: 



for n g]0, 2[. Other values of n could be handled in a similar way, but we focus 
on the range ]0, 2[, which contains n — 1 ("physical" Bures case) and for which 








(9) 
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the results can be expressed without cases discussions. We shall specialize at 
the end our results to n = 1, i.e. for the random Bures states. 

It is not clear to us if Eqn ^ can be derived from a distance on the set of 
density matrices, which could be interesting on its own. For the Bures metric, 
cf. Eqn Q, one can show that the infinitesimal distance is given by (dss)^ = 

Db icr,a + Sa)"^ = | J2i j ^ ^ ^ where the \i are the eigenvalues of a and 
\vi) its eigenvectors [27]. If there exists a distance I?s,„ corresponding to the 
generalized Bures measure, then the associated infinitesimal distance must be 
given by: 

idsB^.r = Db^. (a, a + 5af = \ J] (10) 
1.4 Outline 

Our results on the statistical properties of random Bures states (the case n = 1) 
are summarized in section [2j In section [3] we explain the Coulomb gas and 
saddle point method we use to compute the purity distribution when N oo. 
With this method we derive an integral equation for the eigenvalue density. The 
distribution of the purity for random Hilbert-Schmidt states has been computed 
with the same method [30j. However, in the Hilbert-Schmidt case, this leads to 
integral equations for the eigenvalue density that can be solved using a theorem 
by Tricomi. This method cannot be used in the Bures case as the equations are 
more involved. 

In the distribution of the purity, we find three regimes. In section |4j we 
derive the explicit solution for regime II, i.e. in the case where the density 
has a support of the form [0,5]. One can use for this purpose a theorem of 
Biickner which is the appropriate generalization of the result of Tricomi [SI . In 
section [5] we compute the solution for regime III again using Biickner method. 
This regime is characterized by an isolated eigenvalue (the maximal one) and 
a continuous density for the other eigenvalues. A form of sudden condensation 
happens at the transition between II and III: the maximal eigenvalue jumps 
from a small value to a much larger value, thus detaching from the sea of the 
other eigenvalues. In section [6] we derive the explicit solution for regime I, 
characterized by a density with support [a, b] away from (the most involved 
case). For this purpose, we use a technique developed in the context of the 0{n) 
model by Eynard and Kristjansen |13| . and reformulated in a more algebraic 
geometric way recently by one of the authors and Eynard [SJ 0] . This solution 
involves elliptic functions and Jacobi theta functions. We also analyze the limit 
of low purity (completely mixed) states, and the transition between regimes I 
and II. Some intermediate steps of computations are included in Appendix. 
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2 Statistical properties of random Bures states 

We recall that, for random Bures states, the eigenvalues of the density matrix 
a are distributed according to the law [19l |35] : 



Pb(Ai, ...Xn) = Cn ^ n 



1 li \ i<k ■' 



n,<fc(A,-Afc) 



2 



(11) 



where C^r and D^v are normalisation constants and where the delta function 
enforces the condition = Trcr = 1. 

2.1 Average density and purity 

Sommers and Zyczkowski have computed in [36 the average density of eigenval- 
ues p{X,N) — J2i=i{^ (A ^ Ai)). The quantity p(A,iV)dA is the probability 
to find an eigenvalue between A and A + dA. For a random Bures state, i.e. 
when the eigenvalues Xi are distributed according to the joint distribution in 



Eqn (11), the density of states is given for large N by p{X,N) — Np*{XN) 
where the rescaled density p* has a finite support over ]0, b] with b = 3\/3 and 
is given by 



= ^ h + - --\-.-^] (12) 




The scaling of the density was expected. A typical eigenvalue indeed scales as 
Atyp ~ 1/N for large N because of the constraint J^iLi Aj — 1- 

As we have seen in the introduction, another quantity of interest (and the 
one we focus on in this paper) is the purity defined as: 



N 



Tr[a^]=Y^Xj 



(13) 



E2 measures how mixed (or how pure) a state described by a density matrix 
a = J2i=i Xi\vi){vi\ is. Let us consider two limiting cases: 

• When only one eigenvalue, say A^q is nonzero and thus equal to one, the purity 
is maximal T,2 — 1- In this case, the quantum state is pure, a is the projector 
on state |wio): cr — \vig){vig\. This means that there is no interaction (thus no 
entanglement) with the environment. 

• When all eigenvalues are equal, i.e. Xj — 1/N for all j, the purity is minimal 
E2 = 1/A'^. The quantum state is then completely mixed, it is a mixture of all 
possible pure states \vj) with equal probability 1/A^. In this case, entanglement 
with the environment is maximal. 
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The authors of [5S] compute, for random Bures states, the average of the 
purity S2 — Tr [cr^] = J2i=i ^1 finite N, and the average of the generahzed 
purity = X^ili K large N, and compare their results to the Hilbert- 
Schmidt case. They find (E2) ^ 2^ for large N in the Bures case whereas 
(E2) ^ jf for large N in the Hilbert-Schmidt case. Later, Osipov, Sommers 
and Zyczkowski |31j have shown how to compute recursively the moments of 
the purity at finite N. Though effective, this method does not allow to write 
an analytical answer for the distribution of purity for random Bures states. 

For Hilbert-Schmidt random states, the distribution of (and thus of the 
Renyi entropy Sq — In Sq) was computed for large N and all q > 1 in [^150] . 
It shows interesting phase transitions, that we find again in the distribution of 
purity for Bures random states. Our method could also be applied on to compute 
the distribution of at large N for Bures random states (it amounts to take 
a potential of the form V{x) = + tiX + tqx'^ below), and we expect similar 
phase transitions although the analytical expression of the thresholds will be 
different. 

2.2 Our results: purity distribution for random Bures 
states 

Since a typical eigenvalue is of order (because X^ili ^« — 1); '^'^ ^il^ t)^ 
interested in the behaviour of V (S2 = for s fixed and iV — > 00. To 

compute this distribution, we use a mapping to random matrix theory and a 
Coulomb gas method as explained in sectionjS) Similarly to the Hilbert-Schmidt 
case [35], we find three regimes for the distribution of the purity for random 
Bures states. These regimes follow from two phase transitions in the associated 
Coulomb gas (cf. section [s]) . We show that the distribution of the purity is 
given for large N by: 







exp{- 




for 


1 < S < .Sl 






exp{- 




for 


Sl < S < S2 






. exp{- 




for 


S > S2 



where the critical points si and S2 are given by: 

315 5 
'^= 256' ''^2 

thus S2 I when N ^ 00. The exact mathematical meaning of the sym- 
bol « is a logarithmic equivalent whereas ~ denotes an equivalent: for exam- 
ple r (E2 = jf) ~ exp{-7V2^/(s)} for 1 < s < Sl means InV (E2 = ^) 
—N^^i{s) for large N. Let us define the rate function $ by: 

InV (t.2 = ^) ^ -iV2$(s) as iV ^ 00 (16) 
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$ describes the large deviations of the purity. We show in this paper that 
$(s) = ^i{s) if s < si, <I>(s) = <I>//(s) if si < s < S2 and $(s) = ^ni{s)/\/N if 
s > S2, where ^ijijii do not depend on iV when — > oo. Fig. [l] shows a plot 
of the rate function $(s) (we plot here our analytical results), in particular one 
can see the three regimes. In the three subsections below, we give a summary 
of our results for each of the three regimes, in particular we give the expression 
that we obtain respectively for $//, $/// and $/. 



2.2.1 Regime II 



Regime II corresponds to si < s < S2 with Si = |^ and S2 = | 1 + (277) + 
for large N. It describes the central region, around the mean value of the purity 
s — s — ^. In this regime, the density of eigenvalues is a continuous density 



with support ]0, 6] that vanishes at b but diverges at the origin, see Eqn (22) 



We compute explicitly the rate function (cf. Eqn (82 1): 




(17) 



b — b{s) is the upper bound of the density support. 



Regime II is well defined only for 



315 



< s < 



45 
16 ■ 



Indeed b{s) is real (and 



256' 



positive) only for s < yg, and the density pc becomes negative for s < 
which is non-physical. 

For s < Si with si — |||, the density has a support over [a,b] with a > 0, 
this is regime I. For s > S2 (s2 = 5/2 when N — > 00), we find that there is 
actually another solution that becomes more stable: this is a solution with one 
eigenvalue much larger than the other, described by one single eigenvalue plus 
a continuous density, this is regime III. Therefore regime II is only valid for 

Si < S < S2- 



As the distribution of the purity is highly peaked for large A'', the mean 
value of E2 = s/N is also the most probable value. Thus (S2) ~ for large A'' 



where s minimizes ^//(s), thus 
Sommers and Zyczkowski (36J |^ 



0, cf Fi 



*0 



We recover the result of 





s 


5 


5 


(S2)' 




- — as iV - 


-> 00, le s = - 


^ N ^ 


2N 


2 



(18) 



Around its minimum s — 5/2, the rate function $7/ has a quadratic behaviour 
^//(s) ass^'l, which means that the distribution of the purity 



^Thcy indeed showed that (S2) 



5iV + l 
2JV(JV^+2j 



for large N. 



10 



1.0 - */(s) 




'in lis) 



0.01 













1 */i/(s) 



















25 
S 



Figure 1: Distribution of the purity S2 for random Bures states (n = 1, 
/X = 2/3) for N ~ 5000: (a) Upper figure: plot of the rate function $(s) 



N 



^lnP(E2 = jj,Ny The three regimes are shown. Regime I goes from 



s = 1 to s = si = ^ w 1.23. Regime II goes from s = si to s = S2 ~ 2.70 for 
N = 5000. It includes the mean value s — s — 5/2. Regime III is valid 
for s > S2- The scaling with N (for large N) is different for regime III: 
4>(s) = ^jjj{s)/VN. This is the reason why <I> seems very close to zero within 
regime III. (b) lower figure: zoom of the previous figure around s = S2. At 
the transition I-II, ie at s = S2, we clearly see that the first derivative of $(s) 
is discontinuous. This figure also shows that there is a small range where the 
solutions of regime II and III coexist (solid and dashed lines). The one which 
dominates (solid line) at a given value of s is the one with lowest rate function 
(lowest energy). 
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has a Gaussian behaviour around its average: 



S2 = -j « exp |-— - 2 j I for s close to ^ (19) 

Again, as the distribution is highly peaked for large N, the variance of E2 is 
given to leading order in N by the variance of the Gaussian describing the 
neighbourhood of the mean value: 



2^5 

VarE2-(S2)-(S2)2^^ as iV 



(20) 



Our result is in agreement with [31] . where the authors find for the first mo- 
ment^^i ^ (E^) = ^A^l+l^ and ^2 ^ (Si) = ,^^^.l%;'r,''Zt^l,y With 
their results, one obtains for large N: 

2_ 25iV6 _ 71^4 ^ 7Q^2 _ 24 25 
Var U-^^2-^i^- ^ 2)2(7V2 + 4){N^ + 6) ^ Im ^^^> 

Since the eigenvalues are typically of order 1/iV for large N, the density of 
eigenvalues of a (for fixed purity s) is expected to be of the form pnW = 
Npc(XN) where pc is the rescaled density independent of N for large N. In 
regime II the rescaled density pc has a finite support ]0,b]. It is explicitly given 
by: 



(22) 



with a — (45-\/3 — 166) and (3 = ^{b — 3\/3) where the upper bound of 

the density support is given by 6 = b{s) — |-\/3 ^1 — y^l — . The density 
vanishes at b but diverges at the origin: 

pAx) ~ as a; — 0, PcM ^ C2 Vb — x as x — ?> 5 (23) 

where ci = 2^b'^^^\ ^'^'^ = ^ ^^^^^ ■ ^ plo* of the density for s — 2 
(thus within regime II) is shown in Fig. ^ (blue solid line) . 




In this paper the general formula for /ij, is correct but in the application to fc = 2, there 
is a small mistake, one should read instead of ^(^i''^^)^j^2+i)\^''^+^y 
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2.2.2 Regime III 



Regime III corresponds to higher values of the purity, s > S2 where S2 



for large N. Larger values of s correspond to purer states. 

In regime III (for s > S2), there is one eigenvalue, the maximal one Amax 
that is much larger than the other eigenvalues. Typically one has in this regime 
Amax ^ whereas A^ x ^ for i ^ max. The density is thus made of two parts: 
a continuous density with finite support ]0,6] describing the {N — 1) smallest 
eigenvalues plus a delta peak at A = A„iax- 

Very large values of s correspond to purer states. In particular the limit 
iS* = s/N — >■ 1 is reached when the system is in a pure state, i.e. when there 
is no interaction (thus no entanglement) with the environment. This limit cor- 
responds to the case where only one eigenvalue is non-zero, i.e. A^^ = 1 and 
Xj = for J 7^ iQ. Therefore the maximal eigenvalue is expected to become 
larger and larger when s increases, but what is more surprising is that there is 
one point, s = S2 where the maximal eigenvalue suddenly jumps from a value of 
order 1/N to a much larger value of order \/^fN. This transition is very similar 
to the one observed in the Hilbert-Schmidt case between the second and third 
regime, cf. ref. I5D1. 



In regime III, we show that the rate function has a different scaling with N , 
with 



i.e. $(s) 



/N 



- 5/2 _ 



(24) 



The scaling TV'^/^ and the expression of $/// arc actually the same as in the 
Hilbert-Schmidt case [50] (except for the fact that s = 2 instead of 5/2). 

We have already seen that the solution of regime II is well defined for si < 
s < j| with si = III < s (s = 5/2). Regime III is well defined for s > s. On 
the range s < s < || the two solutions thus coexist. For large iV, the stable 
solution is the one with lowest energy, i.e. lowest rate function $, cf. Coulomb 
gas interpretation section [Sj see also Fig. [l] The transition between II and III 
thus occurs at s = S2 such that N'^^n{s2) = N^^'^^in{s2)- We find: 
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5 
2N 



as N ■ 



(25) 



And at this value, 7V^<i>// = 2~^/^5^/^ iV/^. For very large A'^, the transition 
occurs sharply at S2 = s = 5/2 

At the transition point, the rate function $ is continuous but its first deriva- 



tive is discontinuous. We have indeed 

,1/3 



(!) 



2/3 



,1/3 



and 



1 djvr 
^ ds 



24/352/3 



(^) 



as A^ 
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2.2.3 Regime I 

Regime I corresponds to small values of the purity, 1 < s < si where si = |||- 
The limit s — 1, i.e. E2 — l/iV describes a completely mixed state, thus a state 
that is maximally entangled with the environment. In regime I, the density is 
a continuous density with finite support [a, b] with < a < & that vanishes at a 
and b, cf Fig. |2] (red solid line). 

In this regime, the computations arc more involved. However we are able 
to compute explicitly the Stieltjes transform of the density W{z) — da; ^^3^ 
for z € C \ [a,b], (and then we can come back to the density): we get a quite 



complicated expression given by a sum of elliptic functions, see below Eqn (27) 



We compute the solution parametrically: all the parameters of the problem 



including the rescaled purity s are expressed in terms of A: = a/6, cf. Eqn ( 135 1 
with /x = 2/3, n = l. The expression for the rescaled purity s as a function of 
k can in principle be inverted to get the parameter fc as a function of s. Thus 
all the parameters are implicitly functions of s. 

We obtain the following expression for the rate function ^/(s): 

, X io ^1 5 1 ln2 , , 6^ , , , 

^j(s) = -— - — - where to = -hb - t2 hln& + 2/i - h 

i\ ! 22 442 2 

with h= ^^{^~ ^(^)) ' " / (^(^) ~ \ 

(26) 

and where ii, ^2, b are expressed as functions of k in Eqn (1351, and thus 



implicitly functions of s. The Stieltjes transform of the density is given by 

W{z) = ^^+t2Z + Axo{z)+Bxi{z) (27) 



where A — —t2 , B = —ti + i2 7 •\/3 with 7 given as a function of k in Eqn ( 130 1 
with ^ = 2/3, n — 1, and the special functions xo and xi are given in Eqn ( 126 1 
with /i — 2/3, n = 1, they are expressed in terms of elliptic theta functions. 

At the critical point s = si = |||, the rate function $(s) has a weak non- 
analyticity. It is continuous and twice differentiable with continuous first and 
second derivatives, but its third derivative does not exist. We show that: 



^ , , ^ , , 16777216V2 , ,5/2 



(28) 



More precisely $//(si) = $/(si) and for (5 = si — s = ||| — s we get: 

^ , , ^ , . 64 ^ 8192(^2 4194304^3 ^^.^ ^ ^ 

$//(s) - $//(si) = 1 1 h 0((5 ) as (5 ^ 

^ ' ^ ' 27 2025 820125 ^ ' 



^ , , ^ , , 64(5 8192(52 16777216^/2(5^/2 „^ ^ , 

^i{s) ^ $,(.,) = ^ + ^ + 11390625 + ) ^ ^ « 



(29) 
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Figure 2: Plot of the (rescaled) density Pc{x) at fixed purity E2 = s/N for a 
random Bures state (n = 1, /i = 2/3), for s = 2 s = 1.23 and s = 1.11. The 
blue solid line corresponds to s = 2 within regime II: the density has a support 
]0, 6] and diverges at the origin. The red solid line corresponds to s = 1.11 
(or k — 0.2) within regime I: the density has a finite support over [a,b] and 
vanishes at a and b. The orange dashed line is the density at the transition 
point s = Si — HI « 1.23 between I and II. Here the exact expressions of the 
density are used to plot this figure. 



In the limit s — >■ 1+, which corresponds (for a) to completely mixed states, 
we show: 

<I>Ks) = -^ln(s-l) + 0(l) (30) 
The distribution of the purity E2 vanishes when s — > 1 as: 

■P (Sa = |^,7V) oc (s- 1)T^ ass^l+ (31) 

The probability that a is almost completely mixed is very small. Note that this 
s — ^ 1 behaviour of the distribution of S2 for random Bures states is exactly 
the same as the one of the distribution of S2 for random Hilbcrt-Schmidt states 

m- 
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3 Coulomb gas method 



The goal is to compute the probabihty density function of the purity E2 = X^i 
for generaUzed random Bures states. By definition it is given by: 



(32) 



where the distribution 'PB,ni^i, ^n) of the generalized Bures eigenvalues (the 
usual Bures case corresponds to n = 1) is given in Eqn (|9|: 

Ps,„(Ai, ...A^.) = D^s(j2K- l] ^'^f^^.^f^ (33) 
VT / n„fc(A,+Afc)2 

This computation can be done for large N with a Coulomb gas method, that 
we explain in this section. This technique is well-known in of random matrix 
theory and has been used in various contexts such as the distribution of the 
maximal eigenvalue of Gaussian or Wishart random matrices [HI UHl 132] ; non- 
intersecting Brownian interfaces , quantum transport in chaotic cavities [10] , 
phase transition in the restricted trace ensemble [T or the index distribution for 
Gaussian random fields [U [17] and Gaussian random matrices [25^ . 

3.1 Coulomb gas 



The idea is to see the joint distribution 7-'b,„(Ai, ...Xn) in Eqn (331 as a Boltz- 
mann weight: 



VbAM, -Aw) = ^e-^~[{^'>l where 



Zn 

EN[{X^}] = -2 51 In I A, - A,| + ^ ^ In |A, + Xk\ with J] = 1 

j<k j,k i 



(34) 



The \i can thus be seen as the positions of N particles on a line (Ai G IR+) 
with effective energy i?jv[{Ai}] (and partition function Zat). Then 7'b^„(Ai, ...Xn) 
represents the probability to find these particles at positions A^ in the canon- 
ical ensemble at finite temperature (temperature 1 with kg — 1). The en- 
ergy En consists in a two-body logarithmic interaction potential. The part 
— 2 X^-^j, In |Aj — Afel is the Coulomb interaction in two dimensions. It appears 
very often in random matrix theory (RMT). The mapping from random matrix 
eigenvalues to a Coulomb gas problem is well-known in random matrix the- 
ory [H] and has been recently used in a variety of contexts. 

For large N, because of the constraint J2iLi Ai = 1 we typically expect 
Ai X i where x means that A^ is proportional to 1/A^ when iV — >■ 00. Let us 
introduce the rescaled eigenvalues Xi — XiN that are expected to be typically 
independent of N for large N. Therefore the typical scaling for the purity 
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T,2 = X^i is expected to heT,2 ^ jj for large N. Let us thus write S* = ^ in 



Eqn p2h . The equation ( 34 ) can be rewritten as 



where i5Ar[{xi}] = —2 In \xj — Xk \ + — In \xj + Xfej with xi — N 



2 



and with Z^v = N ^ Zj^ the new partition function. 



(35) 



Let us introduce the rescaled density p{x) = (5(a; — x^). It is normalized 

to one by definition {J^ dx p{x) = 1). When N — > oo, we expect the density 
p{x) to have a finite and continuous limit. We will thus first assume that it 
is the case here. It turns out a posteriori that this assumption can hold for a 
large range of values of the purity S = s/N. However we will see that it cannot 
be true for large values of the purity. We will indeed see that for s > S2 for 
some S2 (which will be shown to be 5/2), there is one eigenvalue that becomes 
much larger than the other: the density thus tends to a continuous part plus an 
isolated delta peak. We will call regime III the regime s > S2 (see section [5| . 

Let us first assume that p{x) tends for large to a continuous limit. This 
will actually give the first two regimes (corresponding to s < 52)- In this case, 
we expect from Eqn ( 35 ) that ^ N'^ for large N and more preciseljQ 



/•oo /'C30 n \ 

8n[{x;}]^N^ dx dx'p{x)p{x')[-\n\x-x'\ + -\n\x^~x'\) (36) 
Jo Ja ^ 2 / 



where ~ means "equivalent" for large N. The integrals run from to +00 as 
the eigenvalues of a density matrix are nonnegative, i.e. > for all i. 

The effective energy £jv ^ N"^ is large for large TV, thus the multiple integral 



in Eqn (32) can be computed by a saddle point method. To leading order for 



large N ^ only the minimal energy contributes to the integral. We get 

7^(22 = ^,^)«^e-^'^^[''^l as TV ^00 (37) 
(« means logarithmic equivalent) where pc minimizes the effective energy Eg 



* Let us introduce the rescaled two-points correlation function p^'^\x,x') = 
jv(jv-i) ~ ~ ^j)- P^^' (^'1 (exactly as p{x)) is expected to have a finite and 

continuous limit when N ^ 00. Therefore we get £jv = N{N — 1) f dx f dx'p'^' {x, a;')(ln \x + 
x'\ — 2ln\x — x'\) + Y f dxp{x) In |2a;| thus £jv " N''^ as A'' — >• 00. Also for large N by definition 
of p(^) one has: p^^'>{x,x') ~ p{x)p(x'). 
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that is given by: 



Es[p\ 



dx 



dx p{x)p{x') 



-to 



dxp{x) — 1 





[/ 




Jo 



In I a; 



dxx p(x) — 1 



In Ix - 









[/ 

Jo 



/ dxx"^ p{x) 
Jo 



(38) 

where to, 1,2 are Lagrange multiphers that have been added to take into account 
three constraints in the minimizatiorj^ respectively the normahzation of the 
density J p — I, the unit sum of the eigenvalues 1 — — J dxxp{x) and the 

fixed value of the purity E2 = s/N, i.e. J dxx^ p{x) — s (it replaces the delta 



function in Eqn (32)) 



The effective energy can be formally written as 



Es[p] = / dx dx'p{x)p{x')Uintix,x' 
Jo Jo 

with C/int {x,x ) = —ln\x — x 



dxp{x)V{x) + cq 



(39) 



I + -ln|x + a;'| , V{x) = to + hx + t2- 



and eg ~ — ^ ~ ^ • The charges of the Coulomb gas (that interact with the 
two body-interaction potential Uint as in Eqn (34)) thus live now in an effective 
external potential V{x) = to + tix + ^2^- The shape of this potential depends 
on the Lagrange multipliers that are not fixed yet, and that will be functions of 
the rescaled purity s. Depending on the value of the rescaled purity s, V will 
have different shapes which will lead to phase transitions for the Coulomb gas 
and thus to different regimes for the purity distribution. These phase transitions 
are shown in Fig. [3j 



In Eqn (37), the partition function Zn can be computed by using the same 
Coulomb gas method as above but without the constraint on the purity. Thus 
-Zjv ~ ' where p* minimizes the effective energy without the constraint 

on the purity, E[p] given by: 



E[P] = / 
^0 



dx 



to 



dx' p{x) p{x' 
dx p{x) - 



In |x 



In I 



+ ti 



dx X p{x) 



(40) 



To leading order for large iV, it is not very difficult to see that E[p*] is actually 
equal to i?s[Pc] evaluated at s s where s is the mean value of the rescaled 
puritjj^ (or also its most probable value as the distribution is highly peaked for 
large N), thus E[p*] = Es[pc] and: 



^-N\E,[p^]-E,[p,]) 



as N 



(41) 



^The choice to write 12/ 2 instead of t2 is arbitrary. 
^We will show later that s = 5/2. 
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where « means "logarithmic equivalent" . 



We need to determine the optimal density pc that minimizes [p\ in Eqn ( 38 1 



Sp(x) 



= 0: 



Pc must satisfy 

'pc{x') (in \x - x'l - I In \x + x'\j = to + hx + 



2 / Ax 
Jo 



, X G Supp[pc] 

. ^^^^ 

where Supp[pc] is the support of the density pc, Supp[/9c] C IR+. Differentiating 
with respect to x we get: 



Q X X 



Q X -\- X 



ti + t2X 



G Supp[pc 



(43) 



where j- is the principal value of the integraQ 

We need to find the solution Pc{x) of Eqn ( [43| that is positive on its sup- 
port and satisfies the three constraints dxp{x) = 1, J^dxxp{x) — 1 and 

dx x^ p{x) — s for a given s > 1. 

The support of the density Supp[pc] C IR+ must be finite. Otherwise, by 
taking the limit a; — >■ cx) in Eqn (|43]) we get a contradiction. The most natural 
finite support we can imagine is a finite segment [a, b]. Let us thus assume that 
the density support is of the form Supp[pc] = [a, ^] C IR+. If a > 0, we expect 
by continuity that the density vanishes at the upper and lower bounds of its 
support, i.e. Pc{a) = = Pc{b)- We will see that this is valid for 1 < s < si 
with si = |||, which we will call regime I. Another possibility is that a — 0. In 
this case the density only needs to vanish at b, i.e. Pc{b) = 0. This is actually 
valid for si < s < S2 as we will see (with S2 ^ | for — > oo) and we will call 
this regime regime II. 

Before coming to the explicit derivation of regime I (section [6| andregime II 
(section l4|, let us present the general setting that we will use in both regimes 



to solve Eqn (43). 



3.2 Resolvent and effective energy 

In the next sections (sections [6] and |4]) , we will explain how to find the density 
Pc solution of Eqn (43) in regimes I and II i.e. for a density with support [a, b] 
with < a < &. The method for the two regimes is quite different but uses in 
both cases complex analysis. In both cases it is easier to compute the resolvent 
(that is the Stieltjes transform of the density) instead of the density. From the 
expression of the resolvent, one can then come back to the density. However the 
minimal energy Es[pc] (and thus the rate function <&) can be computed directly 



from the resolvent (we do not really need the density) , cf. Eqn ( 48 ) below 



^Principal value; d."^ = lim.^o [fr' & + d-'^] 
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Figure 3: Figure showing the (rescaled) density Pc{x) and the effective potentiai 
V{x) = to + ti + 12X^/2 at fixed purity Y,2 — s/N in the three regimes. In 
regime I (left), ie for 1 < s < si, the potential has an absolute minimum at 
a positive x, and the density has a finite support over [a, h] with < a < 5 
(the Coulomb charges accumulate close to the minimum of the potential). In 
regime II (middle), ie for si < s < S2, the potential is monotonically increasing 
for a; > 0. The density has a support on ]0, 6] and diverges at the origin (the 
charges accumulate close to the origin). In regime III (right), ie for s > S2, the 
potential is not anymore bounded from below. The maximal eigenvalue becomes 
much larger than the other eigenvalues (one charge detaches) whereas the other 
eigenvalues remain close to each other (described by a continuous density). 



3.2.1 Resolvent 

By definition the resolvent is the Stieltjes transform of the density: 



W{z) = 



dx 



Pc{x) 



eC\[a,6] 



(44) 



where [a, b] is the density support. 



For a continuous density pc with support [a, 6] satisfying Eqn (43 1 and the 



three constraints J Pc — ^, J dxxpc = 1 and J dxx^pc — s, the resolvent has 
the following properties: 

1/ W{z) is analytic on C \ [a, b] with a cut on [a, b]. 

2/ The jump of W{z) when passing through the cut is given by the density: 

pJx) = — (W(x - iO+) - W(x + iQ+)) = —lmW(x + iO+) , x £ [a, b] 
2m TT 

(45) 

3/ The asymptotic behaviour of W{z) for \z\ — oo is the following^ 

W{z) = 1 + 1 + 4+ (-\ as \z\ ^ oo (46) 



z z^ z--' 
4/ W{z) satisfies the following equation 



W{x + iQ+) + W{x-iQ+) + nW{~x)=t2X + ti , xe[a,b] (47) 



'We have indeed W(z) = ^ Jq°° dxpc(x) + dxxpc{x) + -p- Jq°° dxx^pc{x) + ... 
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Eqn (47) corresponds to the equation of the 0{n) model for a potential 
V'{x) = + ti. The usual Bures case corresponds to n — 1. 



3.2.2 Minimal energy 



We want to compute the distribution of the purity for large N or equivalently 
the minimal energy Es[pc] (cf. Eqn ([4l|)). Es[pc] is given by Eqn pS] ) computed 
for p = pc- The optimal density pc satisfies the three constraints f pc — 1, 
J xpc = 1, / x^Pc = s, thus we have -Es[pc] = /o°° Jg°° dx' Pc{x)pc{x'){— In |a; — 
x'\ + ^Injcc + x'l). Using the saddle point equation ( [42| and the constraints 
again, we get 

\esW]^-%A-'^\ (48) 



To obtain the expression of the Lagrange multipliers io,i,2 and thus the energy 
Eg[pr}, the steps are the following: 

• First we must find the solution W{z) of Eqn (471 for fixed (but unknown) ti^2, 



W{z) must be analytic on C \ [a,h] with a cut on [a, 6]. In general, W{z) will 
be of the form: 

W{z) ^ h{z) + W{z) (49) 



where h{z) is a particular solution of Eqn (471 that we can choose as follows 



h{z) 



2V'{z)~nV'{-z) ti 



Uz 



(50) 



and where W{z) is a general solution (with one cut on [a,b\) of the homogeneous 
equation: 



W{x + iQ+) + W{x - iO+) + n W{-x) = for x e [a, b] 



(51) 



• Then we must impose the asymptotic behaviour of W{z) (cf. Eqn (46)). This 
is equivalent to impose the constraints on the density pc (/ Pc — I, j^Pc — 1, 
/ x^pc = s). We thus get the Lagrange multipliers t\ and ti as well as the 
bounds of the density support a and b as functions of the purity s. 

• We must also compute the last Lagrange multiplier t^. We show below that 



it can be expressed as a sum of integrals of W(^z)^ see Eqn (58). 
• Finally, knowing W[z)^ ti, ti and also to as functions of s, we get the expres- 
sion of the energy Es[pc] = ~~ ^ ~ ^ ■ To recover the usual Bures case, we 
impose n = 1. 



Computation of tg 



to can be computed from Eqn ( 42 ) for x 

fe2 



to = -tib' 



t.-+2 



dxpc{x) 



Inlb- 



ln|5 



(52) 



21 



Moreover, for z ^ [a, 5], we have by definition: 



/'OO / \ 1 /'OO 

W{z)^ dx^^ = — / dxp,{x)\n\z-x\ (53) 
Jo ^ - a; dz Jo 

As pc has a support over [a, b] we thus get dxpc{x) In \z—x\ = da;pc(a^) In \z- 
x\ and thus for Z > b: 

\ dxpc{x)\n\b- x\- dxpc{x)\n\Z-x\^- dzW{z) (54) 

J a J a J b 



or equivalently: 



dxpc{x) hi 



dxpc{x) In 



We can now take the hmit Z oo above, we get: 



11 = dxpc{x) In 

J a 

Similarly one can show that: 

12 = dxpc{x) In 

J a 

therefore we get: 



b ~ X 
b 



b + x 
b 



dz[-- W{z) 



dz I W{z) - - ] (55) 



(56) 



dz W{z) 



to = -t^h~t2 — + (2-n)ln& + 2/i - n/2 
with h= ^ , l2 = J dz (^W{z) - i 



(57) 



(58) 



4 Regime II 

In this section, we assume that the density pdx) solution of the saddle point 
equation (43) has a finite support of the form [0,6], i.e. a — 0. The density is 
expected to be continuous, thus we must have Pc(b) = 0. 

As explained in section |3.2l we need to find the resolvent W{z) that must 



satisfy the conditions 1/-4/ of subsection 3.2.1 The solution can be found using 



a method proposed by Bueckner [9] in the case of a polynomial potential V{x) 
(here we have V{x) = + tiX + to)- 

4.1 Resolvent: Bueckner's method ^ 

We recall that the resolvent is defined as the Stieltjes transform of the density. 



cf. Eqn (44 1. As we have seen in the previous section, the resolvent is of 
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the form W{z) = h(z) + W{z) where h{z) is given in Eqn ( [50| ) and where 
W{z) is a general one-cut solution (cut [0,6]) of the homogeneous equation 
([51]): W{x + zO+) + W{x - iO+) + nW{-x) ^ for x G [a, 6]. 

The solutions of the homogeneous equation Eqn (51 1 form a module of di- 
mension 2 over the ring of polynomials of with real coefficients. In general, 
we have: 

W{z)^A{z^)^^{z)+B{z^)^^{z) (59) 

where A and B are polynomials with real coefficients and where {(/)o,(/'i} is a 
basis of the module. Let us make a change of variables by defining w — w{z) 
and ji (instead of n) 



sm w 



and n = — 2 cos (tt/x) /xgJO, 1[ 



The basis functions are then explicitly given by: 

DS [fi {W + f )] 



00 (2) 



COS 



and 01 (z) = 



sin [fi [w + 1)] 



tan w 



(60) 



(61) 



When \z\ — > 00 we have W{z) ^ i, thus W{z) tends to zero. The polynomials 
A and B are uniquely determined by the condition W{z) = W{z) + h{z) — >■ as 
l^l — >■ 00, where h{z) is the polynomial defined in Eqn ([50]). To determine A and 
B for a given polynomial potential V'{x), we thus need to know the asymptotic 
expansion of 4io{z) and (t)i{z) at 00. One can show that for \z\ — ^ 00: 



00 (-z) = 1 ^ /itan 



/iTT 



01 W 



sm 



b _ 
z 

flCOS 



2 z2 

/iTT 



6 



tan 



/iTT 



r3 



o 



1 



1 



M (m' + 2) 



■ cos 



fin 



r2 



2 

t4 + 2fl^ 



sm 
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b 

z 

fin 



53 

r3 



o 



(62) 



• When |z| — >■ 00, we must have W{z) 
00,1 in Eqn (62), we get A and i3: 



A 



bfit2 



0. Using the asymptotic expansion of 
2t2b 



B 



(2 - n)^/2l^ 



(63) 



We thus have an explicit expression for the resolvent W{z) — h(z) + ^0o(2;) 
B0i(z). 



^In the computations we use very often the following relations obtained from n 
-2cos(TTfi): 

/7r/i\ \/2 — n . f''^ti\ \/2 -I- rt 



tan 
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• We now impose the constraints on pc or equivalently we impose that W{z) = 
^ + ^ + -p- + .... We thus obtain three equations that determine ii, ^2 and b. 
After some manipulations, we get: 



2 (-6 - 3n + 26\/4-7i2//) 



6 (2V4 -n^- 6^(2 - n)) 

53/1(1-^2) 



, /2 + n 



(64) 



12 2fi y 2-n 

The upper bound of the density support is thus exphcitly given by: 



b = b{s) 



3 /2 + n /_ 4 / 2-n\ /i- 



2 - n 



1- Wl 



3 V2 + n/ (1 -M^) 



(65) 



Replacing in Eqn (63) ti and t2 by their expression Eqn (64), we get: 

l2-n 2(2 + ^^2) 



B = 



6(1 + m') 
^2(1-^2)52 V2 + n^(l-M2)fo 
12 24 



(66) 



VT+^(l-/i2)6 V2^M(1-/Lt2)62 



To summarize, the resolvent W{z) is explicitly given in regime II by: 



W{z) + ^ + AM^) + i?0i(z) 

2 + n 2 — n 

, , cos [u (w + f)] , , sin [u (w + f )] 

with </>o z = and 0i z = ^ 

cos [i^J tanu; 

z 1 

where - = —, and n = — 2cos(7r/i) 

sin w 



(67) 



The parameters ti, ^ and i? are functions of = b{s) given in Eqn (64 1 and 



(66), where the upper bound of the density support b{s) is given by Eqn (65) 
For the usual Bures case n = 1, and thus /i = | we get: 

351 22^3 „ 324 , 36^3 

D = 



10&2 hb ' 562 5^ 

_ 3 (-27 + 4%/36) _ 54 (3^3 -b) 

2P ' 5P 

and for the upper bound of the density support: 



6 = 6(.) = ^V3 1-^1-f (69) 
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4.2 Optimal density 

The optimal density Pc{x) is obtained by taking the imaginary part of W (cf. 



subsection 3.2.1 ) 



Pc{x) = lmW{x + iO+) for xe[0,b] 



(70) 



where z = x + iO^ with x e [0, b] corresponds to w = ^ — irj with 77 > 0. More 
precisely, we have - = sinw = cosh 77. Therefore we get: 



X \ X 



=±. = ^±J^_1 for ^^^^_^e[0,&] 
"•^ sm w cosh 77 



(71) 



Using the explicit expresssion of W{z) obtained in the previous subsection, 
we get: 



thus 



-(M+l)r; _ p{p+l)v 

cosh 77 



where cosh 77 = — (72) 

X 



Pc{x) 



&2 



- 1 



/3 x 

TT b 



X V X 
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^-1 



a; V 

i+M / 



X \ X 



62 



(73) 



^2^(2 + fi) 3\/2T^(l + fi) 



/3 = 



(1 - fi)fj.b 
3y/2 - n 



(1-^)^252 
6^/2+^ 



(1-7*2)5 /i (1-7*2) 52 
For the usual Bures case fi — 2/3, n — 1: 



(74) 



3 27 
"= 452(45^- le^') , ^=^(&-3^/3) 



562 



(75) 



.here b = b{s) = |V3 (l - ^1 - ^ 

A plot of the density is shown in Fig. [2] and |3] 



4.3 Minimal energy 

To compute the minimal energy Eg [pd 
to first determine tg. We use the expression in Eqn (|58|) 



f - 1^ - cf. Eqn (|48|, we need 



to = -ii6 - ^^^2 ^ (-2 _ „) ^ ^ 2/1 - nl2 



(76) 
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Using the explicit expresion of W{z) given in Eqn (67), we can compute exactly 
the integrals Ii and I2 defined in Eqn (58), we get: 



2/1 - nh 



3 (2(2 + n) (1 - m') + \/4^M (m' - 3)) 

+ (2 - n)-tE + (2 - 7i)V' ' ^ " ^ 
, 3(n-2) 

l-/i2 



2 

4 vr 
^ ^ 2 



(77) 



1-- ) (-5 + In 4) 



where ipi^) — ^ lnr(a;) is the digamma function and the Euler constant. 
Thus we obtain the exact expression of the minimal energy: 



1 A^ u , 3(2 + n) 5^/4 - 



V4" 



(78) 



.(i-l){....(i^).ln2}.^(2-.) 

Eg [pc] is minimal for h ~ h = ^ \J^^^ which corresponds to s 
As the distribution of the purity is highly peaked for large iV, cf. Eqn (37), the 



2(2+rt)(l-Ai^) 



most probable value s of the rescaled purity is also its mean value. Therefore 
the mean value of the purity for generalized Bures states is given by 



- 2(2 + n) (1 - A^^) 

— as A* — > 00 with s = — 

N 3(2 -n)fi^ 



(79) 



For the usual Bures case n = 1, we recover s — 5/2. 

We can also compute the normalisation Zj\i = E^lp^] (partition func- 



tion), cf. Eqn (41) 



4 2 + n 



fj,\ 2 — n 
(80) ' 

thus the rate function in regime II $7/(s) = Es[pc] — Es[pc] is explicitly given 
by: 



$//(s) 



n \, , 3(2 + n) 5\/4 - (2 - n) 

2-T-^+W^ 25^ + ^ 



7 + 21n 



where b = Ms) — —\—^ — ( 1 
^ ' fi \ 2-n \ 



L 4/2-n\ ; 
3 ^2 + n; ■ 



4(2 + n) 
li^{2-n) 



(81) 
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For n — 1, Eqn (81 ) specializes to 




(82) 



Its minimum is reached at b ~ b ~ S-s/S, i-e. at s = s 



From Eqn (65 1, we see that b{s) is real (and positive) only for s < sq = 
(i^n) ^^2* (i-e.s < 11 for n = 1). Moreover pc becomes negative when 



a + 2/3 > 0, i.e.when 6 < 6i = f 



2+n (1-m) 
2-n 2-/i 



(for n ~ 1, bi = 9\/3/8), which is 



non-physical. Thus the regime can be valid only for b > bi, i.e. for s > Si where 

(83) 





f2 + n\ 


(1-m)'(1 + m)(3-m) 


V2-nj 


^^2(2-^)2 



thus Si = HI for n = 1. The range 1 < s < si on the left of regime II is 
described by regime I, with a density pc with finite support [a, b] with < a < 6 
as we will see. 

It is useful to compute what happens close to s = in order to compare 
later with regime I (for s = s^). In particular, we get: 



7b + V' 



1-p 



In 6a 



/2 + n {I- p) 
2 — n p{2 — p) 



(-3 + lO^i- ll^i^ 



(84) 



12(1-^^)2 

We may also compute the expansion of ^//(s) when 

^ , ^ 2-n (-9 + Up - p^ , 
^//(s) = < ..... + In 



-9 + 14/i - ^2 
2 1^ 24(1 -/.t)2 
{2 - n){3 - p){l + p) (si-s 



2 2 



12(1-^)2 



3 1 - 
{2-n){?.-pf 

Ki-pY 



(85) 

We will see also that, although regime II exists for s in the range [si,So] that 
contains in particular the mean s, there exists also another solution for s > s 
that becomes more stable when s > S2 for some S2 G [s, sq [ (its energy is lower 
energy than the present one). We will show that S2 ^ s as — > oo. This second 
solution is characterized by the fact that one eigenvalue is much larger than the 
other, the density is then described by one single eigenvalue plus a continuous 
part, this is regime III. 
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5 Regime III 



For large s, there is no continuous density pc solution of Eqn (43). In the 
limit s ^ i.e. 5 — > 1, the quantum state becomes pure, which means that 
only one eigenvalue of the density matrix a is nonzero. For s large enough, 
the maximal eigenvalue must detach from the sea of the other eigenvalues to 
become larger and larger, thus it has to be taken into account separately from 
the other eigenvalues. We describe it by Amax = T and the other eigenvalues 
by a continuous density /9jv-i(A) = ^i^max S — K) with support ]0, C] 
(similarly to regime II), cf Fig. [sj 

Exactly as for the third regime in the distribution of purity for random 
Hilbert-Schmidt states [5D] , one can actually show that all the eigenvalues apart 
from the maximal one remain of order l/N as in regimes I and II, whereas the 
maximal eigenvalue Amax is in regime III of order 1/\/N -thus much larger than 
the other eigenvalues for large N . Here, for simplicity, we will not show that it 
is the case (but the proof can be done exactly as in the Hilbert-Schmidt case 
|30j ) ■ we will just assume this behaviour and see that it is coherent. Assuming 
these scalings, the density of the {N — 1) smallest eigenvalues can be replaced 
for large N by the rescaled density p{x) — p^-i (;^) which is expected to 
have a continuous and finite limit when N ^ oo (with support [0, 6] as in regime 
II), whereas A^ax 

This change of behaviour of the maximal eigenvalue between regime II and 
III from 1 /N to 1/ \/N can be understood as a Bose-Einstein type (condensa- 
tion) transition. The purity E2 = J2i ^1 fixed to a value S — s/N. When the 
eigenvalues are very close to each other, they can be described by a continuous 
density p{x), which corresponds actually to regime II, but this regime is valid 
only for s smaller than a certain value sq (cf. previous section), which means 
that the purity S2, given in the continuous case hy J dxx'^ p{x), can not ex- 
ceed the value sq/N. For larger s, there must be (at least jj^one eigenvalue that 
detaches from the other eigenvalues and becomes larger. One can then write the 
purity as = S2 = ;^ / dxp{x)x^ + A^g^^ to leading order for large N. Thus 
-^max ^ ;^ ^ 7^ + -^max' which implies that A^^^^ has to become of same order 
as the purity, thus Amax has to be of order l/ViV. 



5.1 New saddle point 

We assume here as explained above that the eigenvalues are described by an 
isolated eigenvalue Amax = T x 1/\/N plus a continuous density pN~i{X) = 
Iv^ X^i^-^max (A — Ai), which can be written for large N as pN-i{X) — Np (XN) 
where p is expected to have a finite continuous limit as iV — >■ cxi, with a finite 
support [0, b]. Because of the presence of an isolated eigenvalue, the saddle point 
has to be modified. The distribution of the purity E2 — J2i ^1 '^^'^ be written 

^"One can actually show that a case where two or more eigenvalues detach is less favorable, 
cf. [30]. 
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as: 



1 



dt / Vpe 



-N'^Ss [p,T] 



(86) 



where the effective energy is now a function of T and p and is given for large TV 
by: 

pb ph ft \ 

8s[p,T]^ I / dxdx'p{x)p{x') [~hi\x-x'\ + -\n\x + x'\) 
Jo Jq ^ 2 / 



dxp{x) — 1 



/ da; xp{x) + T — 1 
Jo 



(87) 



dx x"^ pix) + N ~ s 



O 



InN 



where to, ti and t2 are Lagrange muhipHers (as in regime II) enforcing the 
three constraints resp. J p = 1 (normahsation of the density of eigenvalues), 
Jxp + T =1 {i.e.J2^ K = 1) and / x^p + T^N = s (i.e.fixed purity E2 = s/N). 

Exactly as in regime II, the functional integral in Eqn (86 1 can be com- 
puted for large iV by a saddle point method -the minimal energy dominates the 
integral. With the correct normalisation, we get: 



,-A'"(£,[pe,Te]-£s[p,,T<,]) 



(88) 



5p 



where {pc,Tc) minimizes £s[PiT], thus 

Ms 

OT lp„T< 

h + t2NTc = 



= 



The saddle point equation = 



reads: 



(89) 



The saddle point equation = 



SEjL 
5p 



is an integral equation on p: 



f dx'pc{x') (in \x -x'\-- In \x + x'\) ^ to + tix + —x^ , a; e [0, b] (90) 
2 / 2 



After differentiating with respect to x, wc find: 



da; 



,Pc{x') 



dx 



iPcjx') 
X + x' 



ti + t2X 



, X 



e [0,6] 



(91) 



Eqn (91) is exactly the same as in regime II, except for the fact that the con- 
straints on pc are slightly modified: here we have J pc — 1, f PcX = 1 — Tc, 
/ pcX^ — s — T^N whereas the Tc were negligible in regime III. 
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5.2 Solution for regime III 

We can use the same Bueckner's method as in regime II to compute the resolvent 
W{z) = J dx^j^ where pc is solution of Eqn (91 1, with the constraints J pc = 1, 



J PcX^l- Tc, J Pc 



T^N. W{z) must be solution of 



W{x + iO+) + W{x -iQ+) + n W{-x) = t2X + h for x € [0, b] 

1 1-T s~ T^N , , (92) 

- H ^ 1 5^ h ... as z oo 

z z-' Z'^ 

The asymptotic condition is equivalent to the three conditions on the density 
Pc (cf. above). The solution of Eqn (92) is as in regime II of the form: 



with W{z) 



W{z) ^ h{z) + A(f>o{z) + B M^) 



(93) 



where h{z) — + is a particular solution, cf. Eqn (50 1, and where 4>q 
and form a basis of solutions of the homogeneous equation, cf. Eqn (67): 



cos 



[^(-+111 and M^) = ^^^[^(- + ^)] 



cos 



tan w 



where 



1 



(94) 



sm w 



and n = — 2cos(7r^) 



When \z\ 
we get: 



cx), we must have W{z) — >■ 0, thus exactly as regime II, cf. Eqn ( 63 ) 
ti bpt2 ^ 2t2b 



2 + n 



B = 



(2-n)V2T^ 



(95) 



We now impose the constraints on pc or equivalently we impose that W{z) 
1 , 1 , s 



when \z\ oo. This step differs from regime II. We get: 
2(2 + n) 



<2 = 



-3 + 3T, + 26,/|^^ 



6(2-n) (2(1-T,) 



(96) 



(97) 



12 (s - iVT^) ^ = 65 



"(l-r,)(l-^2)-6V(l-A^^) (98) 

BE, 



2-n 

Eqn (98) together with the other saddle point equation = 0, ie ti+t2NT^ 



0, give in the large N limit, for > l/N 

I s — s 



Tc 



s 

'3 



Q - 



O 



Ar3/2 



with b = 



N Vs-s 
2 l2 + n 



O 



(99) 



/i V 2 — n 
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where s = | ^§3^^ J2 ^ is the mean value of the rescaled purity, cf. Eqn (79 1 

(for n = 1, i.e. the usual Bures case, we have s = 5/2) and where b is the value 
of the upper bound of the density support at s = s in regime II. Then, using 



Eqn (96 1 and (97) we find the expression of the Lagrange multipliers for large 
N: 



to = 



1 {2-n) 



O 



(100) 
(101) 



The computation of the minimal energy £s[pc,Tc\ is done along the same 



lines as for regime II. It is given by (cf. Eqn ( |87| ) and (90 1): 

£s[Pc,Tc] 



where to = -tib - + (2 - n) Inb + 2h - nh 



(102) 



with h XT {-z - ^(^)) and h = dz {W{z) ~ \). The integrals h 
and I2 can be computed after the change of variable z w where | = • 
After some simplifications, we get when N ^ 00: 



£s[Pc,Tc] = Es[pc] 



(2-n)V^ 



2VN 



(103) 



where Es[pc] is the value of the energy at s = s in regime II, cf. Eqn (80) 

1 — /i 



7b + 



In 



4 2 + n 



Thus £s[pci Tc] = Eg[pc\ and we get for regime III as — ^ cxd: 



p \ 2 — n 
(104)' 



-N2<S>iji(s) 



with $7/7(5) 



(2- 



(105) 



The solution of regime III exists for s > s, cf. Eqn (105). For s < s < Sq 



the two regimes II and III coexist. The regime that is valid for a given s is 
the one that has the minimal energy, cf Fig. [l] The transition between II and 
III thus happens at s = S2 defined by Es^ipc] = £s2[Pc,Tc], i.e. N'^^ii{s2) = 
A^3/2 <^jjj{s2). We find explicitly: 



54/3 



S2 



S + 



with s 



2 (2 



In the limit TV — cxd, s2 tends to the mean value s 



1-^2 



(106) 
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We define the rate function $(s) as 



$7/(s) 



for Sl < S < 32 



*7//(s)/viV for s > S2 



(107) 



At the transition point, the rate function $(s) is continuous but its first deriva- 
tive is discontinuous, cf Fig. [l] 



d$(s) 



ds 
d$(s) 



1_ rf$///(s) 

N 



ds 



(2-n) 



mi st 



ds 



d$77(s) 



ds 



1 _2 



(108) 



The transition between regime II and regime III is thus a first order transition 
(exactly as in the Hilbert-Schmidt case (BO,). 



6 Regime I 

In this section, we assume that the density Pc{x) solution of the saddle point 
equation (43 ) has a finite support of the form [a, b] with < a < &. By continuity, 
the density must vanish at a and b, thus Pc{b) = = Pc{a). 



6.1 Resolvent 

As in regime II, W{z) can always be written: 



W{z) = h{z) + A{z^) xo{z) + B{z^) xi(z) 



(109) 



where h{z) — + is the particular solution given in Eqn (50), and 
where xo ^-iid xi form a basis of the solutions of the homogeneous equation 
W{x + iO+) + W{x ~ iO+) + nW(-a;) = for a; e [a, b\. Now [a, 6] is disjoint 
from [— 6, — o], thus guessing an expression for xo a-nd xi is more difficult. The 
techniques to find a basis (xo , Xi ) were basically developed in [13] , and reformu- 
lated in terms of algebraic geometry in [1]. We first need to introduce elliptic 
functions and theta functions. For a review on those kind of special functions, 
see the book [41^. 



6.1.1 Elliptic parametrization and special functions 

Let us first introduce some useful parameters. We define the modulus fc S [0, 1] 
and the complementary modulus fc' as follows: 



k = 



k' = Vl - fc2 



(110) 
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The complete elliptic integral K and the complete elliptic integral for the com- 
plementary modulus K' are defined by: 



K = 



dt 



^il-t^)il-kH^) 



K' = 



dt 



^{l-t^)il-k'H^) 



(111) 



For k e [0, 1], we have K eR and K' e R. We also define (r, q) and (r', q') such 
that: 

iK ^, _ 1 _ iX' 



K 







q = e"" , 


q' = e*'^^' 



(112) 
(113) 



For k e [0, 1], we have q' e R and g G R. When fc then q' 0. When 
fc -> 1 then g ^ 0. 

We recall n = —2 cos 7r/i, e]0, 1[. Let us introduce the elliptic parametriza- 
tion as in [5]: 



u{z) 



ib 



dz' 



2K' J_, ^(z'2-a2)(z'2-62) 



(114) 



where the integral runs fi'om — a to z in the complex plane following a path 
that has to be specified (u(z) depends on the path). The determination of the 
squareroot k{z) = \/ (z^ — a^)(^z'^ — is choseir^so that it is positive along 
the direct path from —a to points z S] — a, a[. Then, k is an analytic function on 
C \ [—b, —a] U [a, 6], and takes a minus sign across the two cuts [—b, —a] U [a, b]. 

For a path included in H = {z G C, Im z > 0} (for example a path from —a 
to +00 or to — oo of the form z' — y + iO^ with y £ R), we have the remarkable 
values: 



m(— a) = 0, 



u{b) = v!' 



1 



u{go) = u° 



-1 



u[a) = u°- = T 
u{~b) = 



(115) 

However, for a path included in (— BH), for example from —a to +oo (or —a to 
— c») of the form z' = y — iO^ with y G R, we have other values: 



u{—a) — 0, 
u{b) = l+r, 



u{0) = 
u{oo) = 



u{a) — 
u{-b) = 



With a path from —a to —a encircling once the segment [—6, —a] clockwise, we 

1 + i 

2 ' 2 



jet u{-~a) 



1. 



-'^^This can be obtained by writing = —\Jz — ay/z + z — b\/ z + b where the four 

squareroots are chosen with usual determination, i.e. cut on IR_. 
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We may construct a Riemann surface C by gluing two copies of C along the 
segments [—6, —a] and [a, 6], such that the squareroot k{z) defined an analytic 
function on C. The squareroot determination on the second copy is chosen 
opposite to the determination on the first copy, so that k has no discontinuity 
on the [—6, —a] U [a, h] C C. Moreover we have u{+oo) = u{—oo) and actually 
all the points [zj — >■ oo can be identified to one point called oo. Each sheet is 
thus topologically a sphere with two cuts. Therefore the surface obtained by 
gluing the sheets together is topologically equivalent to a torus. This torus has 
two distinguished uncontractible paths. One of them is a path encircling [a, b] 
with positive orientation (or equivalently [—6, —a]). Each time the path from 
—a to z encircles [a,b] (or [—b,—a]), it adds +1 to the value of u{z). Another 
uncontractible path is a path going from —a to +a in one sheet and then from 
+a to —a in the other sheet. This loop adds t + t — 2t to the value of u{z). 
Those two paths form a basis of the homology group of the torus. 

By inverting u(z), we get a function z{u) that is thus doubly periodic with 
periods 1 and 2r (with r G iR'^ here), i.e. z{u + 1) = z{u) and z{u + 2t) = 
z{u). More precisely, one can show that the parametrization has the following 
propertiesp^ 



z{u + 1) — z{u) , z{—u) — z{u) , z{u + r) = —z{u) (116) 



and thus 

z{u) = z{u) , z{t — u) — —z{u) (117) 

where we have defined u as u — 2t — u. 

If we restrict the function u{z) to one upper half plane, one can show that 
u : H — ^ C is a conformal map from BH to the interior of a rectangle with vertices 
— 0, T and — I +T. The inverse function z{u) is thus an analytic function from 
the interior of this rectangle to C. It can be analytically continued on C using 



the symmetries and periodicities of Eqn (116). We thus get a double-periodic 



analytic function u — > z{u) from C to C, namely an "elliptic" function. 

The function z{u) thus belongs to the class of elliptic functions. From the 



definition of u{z) Eqn (114 1, one can show that z{u) is actually closely related 



to the Jacobi sine function sn. Indeed: 

(118) 



z{u) = asnk((f>) where ^ 2) — ' 



z{u) has a pole at — . From its definition in Eqn 114 the expansion 



of z(u) when u — ^ u°° reads: 

z{u) = + xi{u - u°°) + 0((w - u°°f) (119) 



u — 



^"^ The property z(—u) = z(u) comes from the fact that the two sheets represent the complex 
plane but with two opposite determinations of the squareroot. The starting point of the path, 
—a, is the same point on both sheets but the end point 2 £ C can be chosen on one sheet or 
the other. Thus, if for a given path from —a to z we get u = u{z) then the path obtained by 
inverting the two sheets gives u(z) = —u. 
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with: 



ih 

2K' 



bT_ 

2K' 



Xl 



6x-i 



(120) 



6.1.2 Homogeneous equation 



Using the eUiptic parametrization defined in the previous section, we can now 
solve the homogeneous equation 



W{x + iO+) + W{x-iO+)+nWi-x)=0 , xe[a,b] (121) 

For X E [a,b], we have u = u{x + iO+) G r + [— ^,0] for a direct path from 
—a to a; + iO+ (included in H), whereas u{x — ii)^) = 2r — u — 1. Besides, 
u{—x — iO~^) — T ~ u and u{—x + iO+) ~ u ~ t. By abuse of notations, we 
shall write W{u) for W{z{u)). If W{z) is a 1-cut solution (cut on [a,b]) of 
Eqn ([m]), then W{u) + W{2t - u - 1) + nW{T - u) = for u G r + [-5,0]. 
Moreover, W(u + 1) = W{u) and W{—u) = W{u) as z{u + 1) = z{u) and 
z(~u) — z(uYH W{u) is initially defined in the interior of the rectangle of 
vertices —1/2, 1/2, r + 1/2, r — 1/2. We have just seen that its values on the 
boundary of the rectangle are related, and we can use those relations to define 
W{u) as an analytic function for u e C. Then, the relations are equalities 
between analytic functions, so must hold on the whole complex plane: 



Vmgc, n{u + i) = n{u), n{-u)^n{u), 



(122) 



To summarize, to any 1-cut solution W{z) of Eqn 121 we have associated a 
meromorphic function u 1— W{u) defined on the whole complex plane. 
If we focus on the system of equations 



n{u + 1) = il{u), n{u + 2t) + nn{u + r) + Q.{u) = 0, 



(123) 



the space of solutions is of dimension 2 on the field of meromorphic func- 
tions invariant under 1- and 2r-translations. We recall that n = — 2cos(7r/i). 
Then, if we introduce T the operator of r-translation, one can remark that 
ker (T^ -I- riT -f id) = ker (T - e"^id) © ker (T - e-^'^^'id) . A basis of solutions 
is thus given by functions which take a phase e**^^ under r-translation and 
are invariant under 1-translation. We may call "pseudo-elliptic" this kind of 
functions. 

One can show, cf. 13,|, that the unique solution g{u) — Vliu) of Eqn (122 1 



K.{z) 



satisfying the boundary conditions g{u{z)) ^ ^ when z — >■ 00 and g{u{z)) oc 

when z — ;> a, 6, —a, —h where k(z) = ^JJ^^^^^~aF){z^~^^l)^ , is given by a ratio of 
Jacobi theta functions: 



z A (u(z)-§ + f|T) 7?i (-i|r) 
«(z) i),{u{z)-^\r) ^,{^\r) 



(124) 



^•^ However, as W{z) has a cut on [a, b] we have W(x + iO^) ^ W{x — iO'^) thus W{2t — u) ^ 
W{u). 



35 



The Jacobi theta function here is by definition (cf. [H] ) 

1?1 {v\t) = i^(_l)Pe«^^(P+l/2)" + (2p-l)7r« 



(125) 



It satisfies (v + 1|t) = —di {v\t) and di {v + ■ 



-"(^+2^) {v\t). 



g{u{z)) and g{u{—z) ) form a basis of solutions of Eqn 123 By abuse of notations, 
we shall write g{z) for g{u{z)). Notice that g{z) has a zero when z = c, with 
c = a STik{i^,K'). 

g{z) does have a discontinuity on [—5, —a]. Indeed, we have not taken into 
account yet the equation W{u) = W{—u) which mean that W{z) has no dis- 
continuity on [—6, —a]. We can now enforce parity to construct a basis (xoiXi) 
of 1-cut solution of Eqn |121| 



n^{z) f{z) 



and xiiz) 



with k(z) — \/ [zP- — a^)(z^ — b^) 



where / and / are defined by: 



m . + and fiz) = ^ 



where g is the pseudo-elliptic function defined in Eqn (124) 



(126) 



(127) 



6.1.3 Resolvent 

As we have seen above, the resolvent W{z) solution of W{x + iO+) + W{x 
iO+) + nW{~x) = V'{x) = t2X + ti is given by: 



W{z) = h{z) + A{z^) xn{z) + B{z^) xi{z) 



where h{z) is the particular solution h{z) 



2V'(z)-nV'(-z) 



(128) 



2T;r + and 



where xo s-^d Xi ^-re given by Eqn ( 126 ) above. 

A(z'^) and B{z^) are polynomials of that are determined using the asymp- 
totic condition W{z) — >• as |z| — > oo. 

Then we impose the three conditions on the density pc or equivalently 
we impose the asymptotic expansion of W{z) when \z\ — > oo, i.e. W{z) = 
i + -f ^ -I- .... There is also an additional condition coming from the fact 
that W{z) must be analytic at z = ±c (as c ^ [a, h]). These four conditions fix 
the value of ^i, ^2, a and h (as functions of s). 

In order to write down these conditions, we need the asymptotic expansion 
of g{z) when z — )■ oo. It can be computed with a minimal number of new 
notations by first establishing that [5 : 



1 d.g(z) 
7(z) dz 



-^7 - 



zk{z)+ck{c) 



1 Ak{z) 



k{z) 



k{z) dz 



(129) 
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where: 

7= ^(ln^i)'(^-^|r), c = asnk{itiK') (130) 

It is also convenient to introduce A = £ [R!j_. Then, the expansion of g{z) 
takes the form: 

Z Z'^ I z-^ 

+ ^{7'+7(3c^-4(a2 + fo^))-2c^A})l 

1 r ' (131) 

+ ^|9a^ + Sfe"* - 3c* + Qa^b^ - Ga^c^ - &b^c^ ~ 107^(02 + b^) 

+ 6cV+7^-8cV}l+o(l) 
• The resolvent W{z) must be analytic at z = ±c, it imposes: 

A./f±"A(c2) + B(c2)=0 (132) 
V 2 — n 

When |z| — >■ 00, W{z) tends to zero. Therefore we get: 

/2 + n 

A=-t2, B = ~h+t2i\ (133) 

V 2 — n 

When \z\ — )■ 00, VF(z) behaves more precisely as W{z) = 1 + ^ + ^ + .. 
(cf. subsection 3.2. 1[ ). Taking also account of Eqn ( |132[ ), we get the following 



system of four equations that fix the value of the Lagrange multipliers ti^2 and 
the bounds a, b (for a given purity s): 

2 - n V4 - n2 

1 ^2 1 / 2 I 7,2 2 2\ , ^1 



2-n2 ^ ' ' ' yfl-n 



t2 1 



(a* + &* + 7' - - 20^52 + 2c2(a2 + 5^) - 272(a2 + fo^ + , 
(2Ac2 + 7(a2 + 62 ^ 3c2) - ^ ) 



2 - n 

(134) 
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6.1.4 Parametric solution 



We can solve the system of Eqn ( 134 1 parametrically with the variable k = a/h, 
or equivalently with the variable r. 



ti = -(2-n)ai 



SI 



2 + n S2 
2-1163 



to = 



(2-n)^ 51 
2 + n 51 

2 + n 5i52 
2-n ~W 



(135) 



where 



Qfi = A — 7 



= -(l + fc2-c2 + 72)-7A 



52 

53 - ^(7^ + A22) +7 

(^4 



l + fc2 



7A ^2 



3 2 

fc2 A:4- 



24-T2(l + ^') + l2 + 



2 

1 1 
6 



A7(l 



24 12^^ + ^ ^+ 4 



(136) 



Those are expressed in terms of the basic parameters, which are function of k 
(or r) only: 



c = c/b = fc snfc(i/i_ft'') , A = A/6 

7 = 7/&=^(lni?i)'(^^|r) 
We also record the relations [H]: 



(137) 



if ' 



(138) 



Eqn (135) determines the purity s as a function of the modulus k ^ a/b. This 
relation can in principle be inverted to get fc as a function of s, but it seems 
too difhcult to be done in closed form. Finally all the parameters ti, t2, b and 
a — kb are functions of k and thus implicitly of s. 



6.2 Distribution of purity 



From Eqn ( 48 ) , we know that the minimal energy is given by 



Es[pc 



_to _ ii _ h s 
2 2 4 

with to = -tib - t2 — + [n - 2)\nb + 2Ii - nl2 



(139) 
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where Ji = dz - W{z)) and h = J_^dz {W{z) ~ ^ , of. Eqn ([58|. 
Then the distribution of purity is given by: 

^iis)^Es[pc]-Es[pc] (140) 



where Eg[pc] is the constant introduced in Eqn (80). 

We thus need to integrate the resolvent with respect to dz, and the solution 



in the form of Eqn ( 128 ) is not convenient for this task. It is however a good 
idea to use again the elUptic parametrization, and define uj{u) — z'{u)W{u). 
The new function uj satisfies for all it € C 



uj{u) = ^uj(^u), uj{u + 1) = io{u) 
lJ{u) — nuj[u — t) + Lj{u — 2r) = 



(141) 



We justify in Appendix \K\ that, after some algebra, the minimal energy can be 
written: 



Es[pc 



2 4 V 



' 1 '■ 

-'0 



(142) 



where " Jq" is the constant term (i.e. independent of e) in the asymptotic ex- 
pansion of Je when e — >■ 0, the latter being defined by: 



J. 



-f ie 



du(y — ijj{u) — — Ljj{t — u)) 



In 



(¥) 



where X-i is given Eqn (120 1. Another approach to compute the resolvent, pro- 



viding a residue formula for iij{u), was developed in [4]. It allows to perform in 
a systematic way the integration with respect to du on some paths of the com- 



plex plane, thus to carry further the computation of Eqn (142). The resulting 



expression for uj{u) is given in Eqn (148) 



6.2.1 Residue formula for the resolvent 



From now on, we set v 
thaH 



1 — /i, so that n = 2cos7ri^. It was claimed in [4] 



w(uo) = 2 Res Q{uq,u)u){u) 



2 Res 



2V'(x{u))-nV'{-x{u)) 1 \ (144) 



4 - 



2 



:u)) 



x{u) 



-'^^Eqn | |144| can be justified in two steps. First, one notices that Q(uo,u) is constructed 
to be a solution of Eqn l |141| l with respect to the variable uo, which implies that the residue 
expression is also a solution of Eqn | |141| . Second, the behaviour at the poles of the left 
and right hand side of Eqn ( |144| l match, so that the difference is holomorphic on the whole 
complex plane. But when n £j — 2, 2[, the functional relation implies that such functions are 
bounded, and cannot be constant unless 0. Thus, as an application of Liouville theorem, this 
difference vanishes. 
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(145) 



where Res„_>.„oo means the residue when u — >■ This looks like a Cauchy 
residue formula, and the "Cauchy kernel" is given by: 

where Cui''^) = —z'{u + u°°)g{u + can be writterf^ 

It has the following important properties, merely a rewriting of those of g: 

C.(w + 1)-C.H, CAu + T)^e"-Uu), Uu) - - (147) 

u->0 u 

Notice the simple pole of Cu{u) when u — > 0, with residue 1. 

Let us do now the residue computation for V'{x) = t2X + ti. Notice that the 
term contributes for a simple pole with residue —1. Moreover, the term 

involving is a total derivative, thus has no simple pole in its partial fraction 
expansion at infinity. We thus need to evaluate it up to order 0((u — m°°)^^). 
The result is: 

uj{uq) — 2 Res Giua,u) 

2 Res^ ^(^.0, u) ( + - ^ + 0(1)) 

(148) 

where d2G(uo,u) is the partial derivative of G with respect to its second argu- 
ment (here u), and with: 

V2 = T, , vi = 7^ (149) 



where a:_i is defined in Eqn (120) . Then, we want the combination: 

Tl 7X 

-uj{uo) - -iD{t - uo) = -w(uo) + -w(uo - t) 

= -2giuo,u°") + 2vid2Q{uo, u°°) + V2dlg{uo, u°°) 

(150) 

where Q{uq,u) = —Q{uq,u) + ^Giuo ^ t^u). From Eqn (145), and using the 
relations Ciy('u + t) = e*^"^ Ci^(w) and n = 2 cos(7rzy), we find: 

— n, 

G{uo,u) = ~g{uo,u) + -g{uQ - T,u) 

= \ (e-^'^'^Uu^ + u)~ C.("o ~u)+ C.(-uo + u) - e2-'^C.(-uo - u)) 

(151) 



^^We prefer to write instead of C,, to avoid confusion with the usual notation for the 
WeierstraB zeta function. 
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Moreover u°° — thus 2u°° = r mod Z and thus: 



SK, ^°°) = ^ ( - e'^^'^'^Uuo - "°°) - e^-^/^^.^^oo _ j 

dlGiu^, un = ~ ( - e— ''/^Cl^.o - - e-^/2C(^^°° - «o)) 

(152) 

6.2.2 Consequences 

Now, we can perform the integral and find "Jq"- The terms ^?(uo, it°°) 
for k = 1,2 can easily be integrated given the formulas (152 1. We need the 
asymptotic expansion of Ci^(u) when u — > 0, which can be derived from its 
expression in terms of theta functions, or from that of g(z) when z oo: 

(^{u) ^ - + Co + Ciu + o{u) asu->0 (153) 
u 



C.--^. C,^ > f^-c'-y) (154, 



where: 

X-i ' 2x'^i 
We thus find: 

"Jo" = I /l^(^^-iCo) + I ix\C,) + "Jo" (155) 

where remains a piece " Jo" for which we could not find a more explicit form in 
general 

"Jo" = limj^ 2 (e--^/2 / duUu)~e^^''/' / d.. 



2 - n , 
In 



(156) 

Actually, when j/ is rational, "Jo" can be written in closed form, but this 
expression depends on the denominator of v. All the same, g and ^i, can be 
reduced to algebraic functions when v is rational. Those expressions are quite 
complicated to work with, and we have prefered to state results with elliptic 
functions and theta functions, whose origin is more transparent. However, since 
one is particularly interested in n = 1 (corresponding to = 1/3, /x = 2/3), i.e. 
the original Bures measure problem, let us state in that case the result for " Jo" . 
Let us define implicitly fc, K and K' , such that: 



iK iK /i54(0|3r) 



^^w^wr ^-^i^^'* 
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With such definitions, K and K' are the complete elhptic integrals associated 
to the modulus k. Then, we have: 



Cl/3 



2iK' 



1 



1 



%/3- 



1 



sn^{2iK'u) 



snj(2iif' (it + r)) 



snT{2iK' {u + t)) 



The primitive of ^ is |41) : 



dw 



snr(t;) 



= - In 



cn^(u) -dnj(w) 



(u) + dn^(w) 



(158) 



(159) 



We mtroduce the notation Ri = ~ ' ' 4 — '^-ri 

using the properties of Jacobi elliptic functions: 



and after some algebra 



4ln 

n=l 2 



K' 



lni?i 



lni?2 - \/31ni?3 



(160) 



However, we will not use this expression in the remaining of this article. 
6.3 Limit s — t- 1: completely mixed states 

The limit s — > 1 (i.e. S — > 1/A^) corresponds to completely mixed states (mini- 
mal purity). This limit is reached when all eigenvalues of the density matrix a 
are equal. In other words, we have in the limit a ^ b, i.e. A; — > 1 or equivalently 
q = e*'^'^ — )• 0. In this regime, the fast converging series in q defining i?i are 
convenient to find the asymptotic expansion of all interesting quantities. 

6.3.1 Asymptotics of ti,t2,a,b and s 
The basic parameters to compute are: 



ks,Tik{i^iK') = ik 



= -m/ ~M + 2(4 + n)q^ ~ 8(2 + n)q^ + 2(16 + 9n + n^)q* 
V 2 — n L 

- 8(7 + 5n + n2)q^ + 2(2 + n)(24 + 8n + n^)q^ - 8(2 + n)(10 + An + n^)q'^ 
+ 2(128 + 129n + 49n^l0n^ + n'*)q** + o(g^)| 



(161) 
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isnk{ilJ,K') 



U3(0|r)J z9iW2|T)^2(/i/2|r) 



V4" 



' + g^(16 - n^) - 4(12 - + (128 - 13n'^)q-^ 



(162) 



- 12(26 - 3n2)q5 + (704 - 84^^ - n%''' - 4(376 - 4671^ - n%'^ 
+ 3(1024 - 127n2 - 4n'^)g* + o(q^)| 



2 



i(0| 



-^4(0|t)^2 
/3(0h 

= 1 - 8q + 32g2 - 96g3 + 2569"* - 624g5 + U08q^ 
- 3008g^ + 6144q* + o(g^) 
^' = I ^^3(0^) = |{l + 49 + 4g2 + + 89^ + 4q» + o(g*^) j 
(ln7?i)'(://2) 



(163) 



2K' 



I2 + n r 

W — ^ n - 4g + 2(8 - n)q^ - 8(6 - n)q^ + 2(64 - lln - n^)^'' 
V 2 — n L 

- 8(39 - 7?i - n^)q^ + 2(352 - 64?i - lOn^ - n^)q^ - 8(188 - 34n - 6n^ - n^)q'^ 
+ 2(1536 - 279n - Sln^ - lOn^ - n'^)q^ + o(g*^)| 



(164) 



Now, the expansions of the parameters defined in Eqn ( 136 1 read: 



ai = W ^ ^(1 -Aq + 2(8 + n)q^ - 8(6 + n)q^ + 2(64 + lln - n^)q'^ 

V 2 + n L 

- 8(39 + 7n- n^)q^ + (352 + 64n - lOn^ + Ti^)q^ - 8(188 + 34n - 6n^ + n^)q'^ 
+ 2(1536 + 279n - bln^ + IQn^ ~ n'^)q^ + 0(9^) | 

^2 = 4(2 - n)^^q^ - 8q^ + (40 - n)q^ - 8(20 - n)q^ + (556 - 38n + n'^)q^ 

- 8(218 - 18n + n2)g7 + (5056 - 474n - 38^^ + n^)q^ + o{q^)^ 

63 = 4V4-n2|q2 _ 12^3 ^ ^gg _ ^-^^i _ ^^^24 - 3n)q^ + (2348 - 84n + n^)q^ 

- 4(2446 - 112n + 3n'^)q'^ + (36928 - 2006n - 8471^ + n^)q^ + o(g*)| 

54 = 4(2 - n)|q2 - iQq^ + (156 - n)g* - 16(72 - Ji)q^ + (7020 - 150n + n^)q^ 

- 16(2314 - 66n + n^)^^ + (174240 - 6126n - 15071^ + n^)q^ + o(g®)| 

(165) 
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We deduce the expansion of the physical parameters: 



1 8 + 3n 10 + 7n 2 n(n-l) a ,4, 



(166) 



6 = 1 + 4g - 16g^ - 2ng^ + 8(7 - n)q^ + Vlrut + 0(9^) 
a = 1 - 4g + 16g3 - 2ng^ - 8(7 - n)^ + 12ng'' + 0(9^) 
s = 1 + 4g2 - 32g'* + 4(44 - ^n)(f' + o{q^) 

Notice that the position of a < 1 < 6 can be deduced from each other by 
changing the sign of all odd powers in q. This can now be converted to an 
expansion in the variable (s — 1) — >■ 0. We first compute: 

from which we obtain: 

1 3n 9n , , , , , , 
*i = -^-^-^(-l) + o(-l) 

*2 = ^ + - + — s - 1 + o s - 1 

s-1 4 16' (168) 

6=1 + 2\/^ _ J (s - 1)2 + o((s _ 1)2) 

o 

a = 1 - 2x/^ _ ^ (s _ 1)2 + _ 1)2) 
o 

6.3.2 Asymptotics of the purity distribution 

Wc need the asymptotic expansion of Cl/(w) for w pure imaginary between — t/2 
and t/2: 

^i(n-i./2|T) ^^UOk) 

= 7rcotan(7ru) — 7rcotan(7rt'/2) + 47rsin {%{2u — 
+ 47r{ sin (7r(2u - + sin (tt{Au - v)) jg^ (159) 
+ 47r{ sin (7r(2u - 3z/)) + sin ['k{Qu - i>))}q^ 
+ 47r{ sin (7r(2u - Av)) + sin ['k{Au - 2i')) + sin (7r(8u - u)) }q^ 
+ 0{q') 

The 0{q^) is indicated after taking into account the fact that u G [—t/2, t/2] 
and |r| ^00. If u remains bounded, it is in fact a 0{q^). Those expressions 
could be used to obtain an asymptotic expansion of the spectral density. Here, 
we will focus on the purity distribution $/(s), for which we need the expansion 
of "Jo". Since the 0{q^) to integrate in fact reaches the order of magnitude 
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0{q^) from a smaller order of magnitude, and exponentially in the variable 
it remains a 0{q^) after integration over u. 



~ 2 — r? r 

'Jo" = -lnq+ ^— i In 



+ nq^ + (-2+ - +n2)g 



' + o{q')} 
}+oiq') 



(170) 



Notice also that, by expanding Eqn (1691 around u 0, we can identify 
the asymptotic expansion of the constants Cq and Ci (although it could also be 
obtained from the previous expansions and Eqn (154)): 

Co = - 1 - 2(2 - n)q^ - 2(2 - n)(l + n)q^ - 2(2 - n)n\'' + o((z'5)} 

2 

Ci = yl - 1 + 12nq^ - 12{2 - 2n - n^)q^ + Un^q^ + o{q^)^ 

(171) 

Finally, the purity distribution ^/(s) = Es[pc] — Es[pc\ has a small q expansion: 



1 2-nr, /2 12 + n\ /1-M\ 1 1 

l^^'^+4 + — {^-(^V^)+HV^)-^W-2} 



4^ 



8- 7n 



(172) 



/ + o(g^) 



and in terms of the variable s: 



s - 1 



1 2-71 

-4 + — 
-^(-1) 



{ln( 



/i V 2 — n 
104 + 5n 



1-p 



64 



(173) 

The rate function $/(s) has thus a logarithmic (integrable) divergence when s — >■ 
1, with the same leading coefficient —1/2 for all values of n. The distribution 
of the purity thus vanishes rapidly when s — > 1: 



(l]2 = ^,iv) (X(5-1)^ 



as iV — cx) and s 1 



(174) 



The probability that the Bures random state is a highly mixed state, ie s close 
to 1 (it corresponds to the case where the state is highly entangled with its 
environment) is thus very small (it was the same for random Hilbert-Schmidt 
states [50]). 



6.3.3 Remark on separability 

Let us assume that H is the Hilbert space of a bipartite system, i.e. H = 
T-La®T~Lb with N — NaNb (with N the dimension of H, Na the one of Ha and 
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Nb oiHs ) • We recall that a state is said to be separable if it is a convex sum 
of tensor product states, ie if the density matrix is of the form 



(T = ^Picrf (g) erf with Pi > and ^ = 1 (175) 

i i 

and where af G C((H.a), erf € £{'Hb)- otherwise it is entangled. Gurvits and 
Barnum j3] showed that states with purity S2 smaller than jjzri are necessarily 
separable (independently of the bipartition A,B). In our notations, this cor- 
responds to states having s < s' = j^^, which is close to 1 for large N. We 



can deduce from Eqn ( 173 1 a logarithmic equivalent for the probability that a 



random Burcs state has purity less that jqzzi' 

pf-<— UiV-* (176) 

Since the probability that a random Bures states is separable is of order N^^ 1^ 
(Corollary 2 in [12]), we have: 

This means that the criterion of separability based on purity is superexponen- 
tially loose when N is large. 



6.4 Transition I-II 

The transition between regime I and regime II happens at s = si when a — > 0, 
i.e. r or equivalently 1, (( — e*^"^ — e^*'^/'^ ^ 0. It is now convenient 
to perform a modular transformation on the theta functions, to write them as 
fast converging series of g' instead of q: 

^Mr) = (-«t)-i/2 e-'""'/" ^i{v/t\t'), t' = 1/t (178) 

It turns out that an expansion up to 0{{q'y^'^) of all parameters is necessary 
to find the order of the transition between I and II. 



6.4.1 Asymptotics of ti,t2,a,b and s 

{q')^ and Q — {q'Y^^ are both involved in the asymptotics when q' 0. For 
non-rational /i, they are algebraic independent, so all the terms of the form 
{q')^^Q^ are of distinct order of magnitude. For rational yit, one will have to 
several collect of these terms to find the asymptotic expansion at a given order. 
Since we assume 1/2 < /i < 1, and we want expressions up to o((q')^"'"'' = 
{q'Y^Q), we can use (g')^'' S o((g')'^+^) to expand all expressions up to order 
3 in the variable (g')'', ^'^^d then find the general coefficient of the expansion in 



46 



Q. We can also simplify expressions since {q'Y^Q G o{{q'Y^^) and {q'Y^Q^ S 
o{[q'Y~^^) whenever j >2. 



^?2(0|r') ^^{^1T' /2\t') 



X = 



j>0 j>l 

1 \\MO\t')) Mf^r' 



1 + E 2J'Q' + 2(g')'' [1 - 2Q - - E + 2(9')'^ [1 - 5g] + 2(9')''' + o((g')'' 



i>3 



J>2 



f ^3(0|r')-|(l + 4(g'rg + o((g')'+^) 



(179) 



Then, we deduce the expansion for the parameters defined in Eqn ( |136 ): 
ai = 1 - M + 2{qr [1 - 2(1 - m)Q - + 2(9')''^ [l - 4Q] + 2(g')''^ + o((q')'+^) 
62 = + 2(g')^ [ - A* + 2/i(2 - /i)Q - (2 - fi)Q^] + 2f,{q')'^ [-1 + 8Q]- 2fi{q'f'^ 



o{{q') 



A^(2-m)(1-m) 
6 



+ {qT [ - /i' - 2m(2 - /^)(1 - fi)Q + (2 - m)' 



^i'{qr'' [ - 1 + 12Q] - m'(<z')'^ + o{{qr^n 



24 3 



[ - m(1 + m) - 2Ai(l - m)(2 - m)Q + (2 - m)(3 - 



[(1 - m)(1 + A^) + 8(2 + /.2)Q] - Ml m)(1 + m) (^,)3. 

(180) 



We then find that s reaches from below the limit: 

2 + n 3 (1-Ai)^(l + M)(3-/^) 



Sl 



2-n 4 



/i2(2-/x)2 



(181) 
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which gives the value of Si . The hmit of the other physical parameters read: 

- 3 (1 - 

^ (2-n)^ 2 /x^(2-/.)^ 
2 2 + n 9 (1-m)4 (182) 

6* = ./P^3. 



2 — n /x(2 — /x) 



a* = 

and their asymptotic expansion is: 



( V) = (2-M)(l -M)V [^^^ + - - + - ^^(^ - ^)«^] 

/I (-^^2^ 

[(I + + 2/i) + 16(-l - 3m + 



(l-M)4(2-/x) 
4(]^^z)^(lJ^ 
^ (1 - /i)6(2 - /i) ^'^Uyj J 



12(9') 



[m' - 2^(2 - m)Q + (2 - m)'Q'] 



y (2-/i)(i-/i)v 

+ (2-1!)^^-"/^)^ ^ "^^^ " "^^'^^ " "^^^ ~ ^^^^ 

+ 4M1 + m);(6 + 27, + Y-9m3) (^,^3. ^ 



(2-Ai)2(l-/,)6 

b — b* \ 2 fa')'^ 

^~6^) = M(2-M)(l-/x)^ [ ~ '"'^^ + '"^ + ^'"^^ ~ -"^^^ ~ + '"'^^ '"^^^ ~ ^^'^^ 



(2-/^)2(1-^)3 

4M(l + Mf(2 + /x)^ 
(2 - /x)3(l - /x)4 



[ - (1 + fif{2 + 4(2 - /z)(5 + n^)Q] 



(^) = (3-.)(2-.)(l-.)2.(l + .) [ - '^'(^ + + ^^^(2 - '^^^^ - + 



(3-/x)(2-/*)(l-/x)Ml+/") 

-(2-M)2(3-/i)V 

8/i (g')^** 



[(l + M)'+4(2-/x)(l-5M)g] 
(<7')"' + o((g')'+'') 



(3-m)(2-mP(1-/x)3 

4/t(l + /^)4(2 + /^)(2 + 3//,) 

(3-m)(2-/x)3(1-m)4(1 + /i) 

(183) 
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This could in turn be converted to expansions in the variable (s* — s) — ;> 0. We 
just notice here that at leading order: 



2 - n (2 - ^if^i 



2 + n + 

and since q' > 0, it is clear that the regime I must correspond to s < s* 



(184) 



6.4.2 Asymptotics of the purity distribution 

We need the asymptotic expansion of C.^{tw) for w real valued between —1/2 
and 1/2: 



(185) 

By expanding around w — > 0, we also identify the asymptotic expansion of the 
coefficients Co and Ci: 



Co 



itt/t 1 - {q'Y-^^ 

Ci 1/1 



7r2/r2 



{\ A^') + 2(1 - H) Y^^^^, + 2(1 + {q' f^^ + o{{q'Y+n 



(186) 

We also deduce the asymptotic expansion of "Jo" (the derivation is presented 
in detail in Appendix p|): 



■Jo 



2 — n /"Kb 
In 



4-^2 



■ TT 2 — n 



V.(i) - V 



1-M 



(187) 



where V' — (InT)' is the digamma function, and b and X must also be expanded 
according to Eqns (1791 and (183). Notice that the terms of "Jo" involving tt 
and ip are also present in Es[pc] given in Eqn (80). So, they cancel each other 
in $/(«) — Es[pc] — E-s[pc]- We thus find the small q' expansion: 



2 - n r -9 + 14^ - 



24(1 - py 



+ ln 



(2n)(^ ^ p^l + pr _ 2^(9 _ ^ ^ (3 _ ^)2(2 _ ^)q2] 

3(1 — A^) V 2 — fi 
2{2 -n)p{q'Yf' 



[(l + M)'(-l-4/. + Ai'+/.3) 



3(2-m)2(1-A*)''(1 + m) 
+ 4(2 - /i)^(-13 + 5^^ - llAi^ - 5/?)g] 

(2 - n)A*(l + pf{A + 28m + 7a*' - 31^*^ + 5^*^ + 3^*^) ^^ 3^ 



3(2-m)3(1-m)s 



o{{q') 



49 



If we turn this expression into an expansion in the variable (s* — s), we notice 
that up to a global factor, the coefficient of (g')^ and of {q'YQ are precisely 
those found in the expansion of s* — s. This means that $/(s) contains a linear 
term in (s* — s), and no (s* — s)^/^ as we would have naively expected. The 
next correction comes at order (g')^^ ~ (s* — s)^, and going one step further, we 
find a singular term of order {q'Y'^Q = {q'Y^^ ^ (s* — s)^+^/^. More precisely: 

-9 + IV - , 2 - /i^ 



^ (2-n)(3-M)(l + A^)^s*-g ^ ^ (2^n)(3-A^) 



s — s 



12(1-^)2 \ s* ) 8(1-Ai)2 V s* / (188) 

32(2-n)/i2 /(3-/^)(2-/i)(l-/i)2 s*_s^i+i^ 



(l-/x)5(l + A*) V 4a*(1+/z) 

+ o((s*-s)i+i/M) 



In fact, this expansion matches up to 0{{s* — s)^) that of $//(s), see Eqn (85), 
and the transition between regime I and II is characterized by: 

^,U)-^Ms)^ 32(2-n)M^ / (3-M)(2-M)(l-M)^ f:^y+V^ .^39. 

' ^ (i-m)5(i + m)v Ml + ^Ji) s* J ^'^^^ 

So, the purity distribution and its first and second derivative are continuous 
when s — > s^, but its third left derivative blows up. The transition I-II is 
therefore of order 1 + arccos(«/2) ^'^^ n e]0, 2[. 

Acknowledgments 

Part of this work was conducted during the session "Random matrices and 
vicious walkers" at I'Ecole des Houches, to which we are grateful both for at- 
mosphere and funding. G.B. thanks Guillaume Aubrun for discussion and ref- 
erences, and benefited from the support of the ANR project "Grandes Matrices 
Aleatoires" ANR-08-BLAN-0311-01. C.N. thanks Satya N. Majumdar, Gregory 
Schehr and Pierpaolo Vivo for useful discussions and references. 



A General expression for the minimal energy 

We start from the expression: 

i^^.H^-j-f + (-/!+ 2^2- ^hr6+^ + ^) (190) 
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where: 



b 



b 



dz 

oo 

_ b 

lim 

M— H-oo 



lim 



/ , , n , ,\ (2 - n) ln(Af/6) T/(M) -i V{b) 



M 

Mb 



Mi, 



= lim i / du 



(2 - n) ln[a;(u°° + ie) /&] F(x(u°° + ie)) ' 



t2h'^ tib 
4 2" 

(191) 

We have written V{x) = *^ — h tix. Notice that V{x{u°° + ie)) admits an 
expansion in the limit e — > with a constant term due to the x"^, and equal to 
t2XiX-i, where: 

x{u°^+5) = + xi5 + 0{5'^) (192) 
is the asymptotic expansion of x near infinity. Moreover, we know xi from 



Eqn (120). If we introduce 



J.= P dw(-cJH-^cJ(r-u)) -^ln(^) (193) 



and " Jq" the constant term in the asymptotic expansion of when e — )■ 0, we 
can write: 

E.iPc] - - | " - ^(«' + + "^o" (194) 

B Transition I-II: expansion of Jq 

We want to compute in an asymptotic expansion when g' — > 0: 

i7n./2 _ ~'im./2 fl/2 .-1/2 

Je = ^ \e-^^-n / C.(™)d(™) - e-''/2 / C.(™)| 

2 — n f x^i 
In 



2 V 

(195) 

and in particular the constant term " Jq" in its asymptotic expansion when 



e — >■ 0. We have seen in Eqn ( 185 ) that at the transition between regime I and 
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II, 



lim t(^^{tw) — (2i7r) 



,2iTTW 



- 1 



(196) 



and all the other corrections to (^^(tw) are regular when w ^ 0. So, the compu- 



tation of the subleading terms of " Jq" (written in Eqn ( 170 )) is straightforward, 
and we shall concentrate on the leading term, i.e. lims_>.s. " Jq". 

We perform the change of variables y = e^*'^™ and introduce contours as 
in the picture (the quarter of circles C± are both defined with anticlockwise 
orientation). We now have an integral on a portion of the unit circle in the 
2/-plane. Notice that the function C.iy{u{y)) has a cut on the negative real axis, 
and going from above the cut to below the cut along the unit circle correspond 
to send m to m + r, thus for any y £ IR_: 



(197) 



We have: 



y = e 



2inu/T 



plane 
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2 — n , / 6t \ 



|g-W2^ dyy^_ l^'-\^\dyy^_ 1^ 
(f Ayy'^'^-^ r^AM dyy''/^-^ d;/ (-6-^2) |^|^/2-U 



2-n, / 5t 



\2ieK'J 



= -i^^^ - ^{ X + 1- (^) } + 

(198) 

where the o(l) stands for a quantity going to when e — >■ 0. Thus: 



„.=-v/J^J-^{M2) + /d„»l^} (199) 

The last integral is elementary: 

/ dy^ — = - hm { / dyy'^/^-'il-y)-'^- / dy(l-y)-«) 

= - lun / r(./2)r(i-a) _^N 

a-i-i- Vr(z//2 + l-a) 1-a/ 



(200) 

Since r(l - a) = + 7^ + 0(1 - a) and r'(l) = -7^ while r(l) = 1, we 
obtain: 

where V' = (InT)' is the digamma function. 



53 



References 

[1] G. Akemann, G.M. Cicuta, L. Molinari and G. Vernizzi, "Compact support 
probability distributions in random matrix theory", Phys. Rev. E, 59, 1489 
(1999), cond-mat/9809270. 

[2] G. Aubrun and S.J. Szarek, "Tensor products of convex sets and the volume 
of separable states on TV qudits", Phys. Rev. A, 73, 022109 (2006), quant- 
ph/0503221. 

[3] H. Barnum and L. Gurvits, "Largest separable balls around the maximally 
mixed bipartite quantum state", Phys. Rev. A, 66, 062311 (2002), quant- 
ph/0204159. 

[4] G. Borot, PhD Thesis, Chapter V, Universite d'Orsay (2011), math- 
ph/1110.1493. 

[5] G. Borot and B. Eynard, "Enumeration of maps with self avoiding loops 
and the 0{n) model on random lattices of all topologies", J. Stat. Mech.: 
Theor. Exp., POlOlO (2011), math-ph/0910.5896. 

[6] S.L. Braunstein and CM. Caves, "Statistical distance and the geometry of 
quantum states", Phys. Rev. Lett. , 72, 3439-3442, (1994). 

[7] S.L. Braunstein, "Geometry of quantum inference", Phys. Lett. A, 219, 
169-174, (1996). 

[8] A.J. Bray and D.S. Dean, "Statistics of critical points of Gaussian fields 
on large-dimensional spaces", Phys. Rev. Lett., 98, 150201 (2007), cond- 
mat/0611023. 

[9] H.F. Bueckner, Journal Math. Analysis and Applications, 14, 392 (1966). 

[10] D. Bures, "An extension of Kakutani's theorem on infinite product mea- 
sures to the tensor product of semifinite ■u;*-algebras" , Trans. Amer. Math. 
Soc., 135, 199-212, (1969). 

[11] V. Cappellini, H.J. Sommers, and K. Zyczkowski, "Distribution of G con- 
currence of random pure states". Physical Review A, 74, 062322 (2006), 
quant-ph/0605251. 

[12] D. Dean and S.N. Majumdar, "Large deviations of extreme eigenval- 
ues of random matrices". Physical Review Letters, 97, 160201 (2006), 
cond-mat,stat-mech/0609651 ; "Extreme value statistics of eigenvalues of 
Gaussian random matrices". Physical Review E, 77, 041108 (2008), cond- 
mat.stat-mech/0801.1730. 

[13] B. Eynard and C. Kristjansen, "Exact Solution of the 0{n) Model on a 
Random Lattice", Nucl. Phys. B, 455, 577-618, (1995), hep-th/9512052. 



54 



[14] P. Facchi, U. Marzolino, G. Parisi, S. Pascazio, and A. Scardicchio, "Phase 
transitions of bipartite entanglement", Physical Review Letters, 101, 50502 

(2008) , quant-ph/0712.0015. 

[15] P.J. Forrester, Log-Gases and Random Matrices, Princeton University Press 
(2010), http://www.ms.unimelb.edu.au/ matpjf/matpjf.html. 

[16] P. Di Francesco, P. Ginsparg, J. Zinn-Justin, Phys. Rept., 254, 1 (1995), 
hep-th/9306153. 

[17] Y. V. Fyodorov and I. Williams, "Replica symmetry breaking condition 
exposed by random matrix calculation of landscape complexity" , Journal 
of statistical physics, 129, 1081-1116, (2007), cond-matdis-nn/0702601, 

[18] O. Giraud, "Distribution of bipartite entanglement for random pure states" , 
Journal of Physics A: Mathematical and Theoretical, 40, 2793 (2007), 
quant-ph/0611285. 

[19] M.J.W. Hall, "Random quantum correlations and density operator distri- 
butions". Physics Letters A, 242, 123-129 (1998), quant-ph/9802052. 

[20] P. Horodecki, M. Lewenstein, A. Sanpera, K. Zyczkowski, "On the volume 
of the set of mixed entangled states", Phys Rev A, 58, 883 (1998), quant- 
ph/9804024. 

[21] R. Josza, "Fidelity for mixed quantum states" , Journal of Modern optics, 
41, 2315-2323 (1994). 

[22] I.K. Rostov, Mod. Phys. Lett. A, 4, 217 (1989); Phys. Lett. B, 266, 312 
(1991). 

[23] E. Lubkin, "Entropy of an n-system from its correlation with a k-reservoir" , 
Journal of Mathematical Physics, 19, 1028, (1978). 

[24] S.N. Majumdar, O. Bohigas, and A. Lakshminarayan, "Exact minimum 
eigenvalue distribution of an entangled random pure state". Journal of 
Statistical Physics, 131, 33-49, (2008), cond-mat.stat-mech/0711.0677. 

[25] S.N. Majumdar, C. Nadal, A. Scardicchio and P. Vivo, "The index distribu- 
tion of Gaussian random matrices". Physical Review Letters, 103, 220603 

(2009) , cond-mat.stat-mech/0910.0775 ; "How many eigenvalues of a Gaus- 
sian random matrix are positive?". Physical Review E, 83, 041105 (2011), 
cond-mat.stat-mech/1012.1107. 

[26] S.N. Majumdar and M. Vergassola, "Large deviations of the maximum 
eigenvalue for Wishart and Gaussian random matrices". Physical Review 
Letters, 102, 060601 (2009), cond-mat.stat-mech/0811.2290. 

[27] C. Nadal, PhD thesis, Chapitre 8, Universite Paris-Sud XI (2011). 



55 



[28] C. Nadal and S.N. Majumdar, "Nonintersecting Brownian interfaces and 
Wishart random matrices", Physical Review E, 79, 061117 (2009), cond- 
mat.stat-mech/0903.1494. 

[29] C. Nadal, S.N. Majumdar and M. Vergassola, "Phase transitions in the 
distribution of bipartite entanglement of a random pure state" Phys. Rev. 
Lett, 104, 110501 (2010), cond-mat.stat-mech/0911.2844. 

[30] C. Nadal, S.N. Majumdar and M. Vergassola, "Statistical distribution of 
quantum entanglement for a random bipartite state", J. Stat. Phys., 142, 
403 (2011), cond-mat.stat-mech/1006.4091. 

[31] V.Al. Osipov, H.-J. Sommers and K. Zyczkowski, "Random Bures mixed 
states and the distribution of their purity", J. Phys. A: Math. Theor., 43, 
055302 (2010), cond-mat.stat-mech/0909.5094. 

[32] D.N. Page, "Average entropy of a subsystem", Phys. Rev. Lett., 71, 1291- 
1294, (1993), gr-qc/9305007. 

[33] A. De Pasquale, P. Facchi, G. Parisi, S. Pascazio and A. Scardicchio, "Phase 
transitions and metastability in the distribution of the bipartite entangle- 
ment of a large quantum system", Phys. Rev. A, 81, 052324 (2010), quant- 
ph/0911.3888. 

[34] A. Peres, "Separability criterion for density matrices", Phys. Rev. Lett., 
77, 1413 (1996), quant-ph/9604005v2. 

[35] H.-J. Sommers and K. Zyczkowski, "Bures volume of the set of mixed quan- 
tum states", J. Phys. A: Math. Gen., 36, 10083 (2003), quant-ph/0304041. 

[36] H.-J. Sommers and K. Zyczkowski, "Statistical properties of random den- 
sity matrices" , J. Phys. A: Math. Gen., 37, 8457 (2004), quant-ph/0405031. 

[37] S.J. Szarek, "Volume of separable states is super-doubly-exponcntially 
small in the number of qubits", Phys. Rev. A, 72, 032304 (2005). 

[38] V. Vedral and M.B. Plenio, "Entanglement measures and purification pro- 
cedures", Phys. Rev. A, 57, 1619 (1998), quant-ph/9707035. 

[39] P. Vivo, S.N. Majumdar and O. Bohigas, "Large deviations of the maxi- 
mal eigenvalue in Wishart random matrices", Journal of Physics A: Math. 
Theor., 40, 4317 (2007), cond-mat |Stat-mech /070137l| 

[40] P. Vivo, S.N. Majumdar and O. Bohigas, "Distributions of conductance 
and shot noise and associated phase transitions", Phys. Rev. Lett., 101, 
216809 (2008), cond-mat.mes-hall/0809.0999 ; "Probability distributions of 
linear statistics in chaotic cavities and associated phase transitions", Phys. 
Rev. B, 81, 104202 (2010), cond-mat.mes-hall/0909.2974. 

[41] Z.X. Wang and D.R. Guo, Special functions. World scientific (1989). 



56 



[42] D. Ye, On the Bures volvme of separable quantum states, J. Math. Phys. 
50, 083502 (2009), quant-ph/0902.1505. 

[43] K. Zyczkowski and H.-J. Sommers, "Induced measures in the space of mixed 
quantum states", Journal of Physics A: Math. Gen., 34, 7111 (2001), 
quant-ph/0012101. 



57 



